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P   procedure  for  automated  design  of  monolithic  integrated 
circuit  amplifiers  is  described.   This  procedure  yields  a 
circuit  which  satisfies  both  ac  and  dc  design  specifications; 
the  ac  specifications  include  requirements  on  input  and 
output  impedances,  gain,  and  bandwidth,  while  dc  specifications 
include  requirements  on  temperature  stability  and  power 
dissipation.   The  design  problem  is  transformed  into  a  math- 
ematical prograiTiming  problem  by  embodying  the  specifications 
in  a  flexible  least  pth  performance  function.   All  ac  and  dc 
specifications  are  considered  simultaneously  because  the 
contributions  to  the  performance  function  from  all  specifi- 
cations are  si:immed  to  form  a  concise  measure  of  the  network's 
ability  to  exhibit  the  desired  characteristics.   Gradients, 
which  are  determined  by  the  adjoint  network  method,  are  used 
to  adjust  the  network  parameters  and  topology  so  as  to 
minimize  this  performance  function.   Topology  adjustment  is 
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achieved  by  growing  and  removing  circuit  elements,  including 
transistors.   Transistor  growth  and  removal  is  indicated  by 
gradient  expressions  that  determine  where  ac  gain  is  necessary. 
Since  each  transistor  can  be  connected  in  one  of  three  config- 
urations, the  most  appropriate  configuration  must  be  selected 
through  use  of  the  gradients.   A  computer  program,  FROLIC 
(Fabrication  Related  Optimization  of  Linear  Integrated 
Circuits) ,  has  been  written  which  implements  the  design  pro- 
cedure, and  several  examples  are  presented  to  demonstrate  its 
ability. 
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CHAPTER  1 
INTRODUCTION 

In  recent  years,  the  use  of  linear  integrated  circuits 
has  increased  significantly.   Development  of  these  circuits 
is  complicated  by  their  compact  nature  and  the  presence  of 
numerous  nonlinear  elements  such  as  diodes  and  transistors. 
In  the  past,  the  circuit  designer  typically  would  build  a 
breadboard  of  his  circuit  to  test  and  evaluate  the  circuit's 
performance.   Often  individual  elements  would  be  changed  to 
improve  the  performance.   With  an  integrated  circuit,  where 
the  entire  circuit  is  on  a  small  chip,  this  procedure  is  not 
feasible  because  of  the  amount  of  time  and  money  required  to 
build  each  circuit.   Therefore,  a  computerized  procedure 
which  automatically  evaluates  an  integrated  circuit's  per- 
formance and  optimizes  it  without  requiring  the  circuit  to 
be  built  furnishes  a  substantially  less  expensive  design  in 
less  time. 

Background 

Network  analysis,  performance  evaluation,  and  parameter 
adjustment  are  three  basic  functions  of  an  automated  network 
design  algorithm.   Network  analysis  is  used  to  determine  the 
responses  of  a  network  to  a  given  set  of  excitations.   A 
scalar  perforraance  function  is  then  used  to  compare  the  actual 


responses  with  the  desired  responses.   Systematic  adjustment 
of  the  network's  parameters  optimizes  this  performance  function 
and  therefore  obtains  the  network  that  best  satisfies  the 
design  specifications.   Each  of  these  operations  must  be  per- 
formed accurately  and  efficiently. 

Design  is  based  upon  optimization  of  the  scalar  perform- 
ance function.   Two  types  of  optimization  methods  can  be  used 
to  adjust  methodically  the  network's  parameters:   direct 
methods  and  gradient  methods.   Direct  methods  which  optimize 
functions  v/ithout  calculation  of  parameter  gradients  generally 
converge  very  slowly  to  an  optimal  set  of  parameter  values 
when  laore  than  five  parameters  are  considered.   Gradient 
methods,  potentially  faster  than  direct  methods,  require  the 
determination  of  the  gradients  of  the  performance  function 
v;^ith  respect  to  all  parameters.   The  classical  technique  used 
to  determine  gradients  is  differencing  which  requires  pertur- 
bation of  each  parameter,  one  at  a  time,  and  reevaluation  of 
the  network's  behavior  after  each  perturbation.   Because  of 
the  large  number  of  analyses  required  using _ perturbation, 
this  method  is  not  much  faster  than  direct  search  methods. 

Prior  to  1967,  computer-aided  design  algorithms  were 
generally  limited  to  network  analysis ^~'*  because  of  the 
large  amounts  of  computational  time  that  were  required  to 
adjust  network  parameters.   A  visual  display  which  interacted 
with  an  analysis  program  was  often  used  to  allow  a  circuit 
designer  to  adjust  network  parameters  and  determine  the  effect 
of  the  adjustment  on  the  network  responses.   Through  a  trial 
and  error  process,  a  final  netv7ork  v^as  obtained. 


In  1967,  Rohrer^  proposed  a  method  for  the  efficient 
determination  of  partial  derivatives  of  a  scalar  performance 
function  with  respect  to  the  network  parameters.   This  work 
v/as  oriented  toward  lumped,  linear,  time-invariant  network 
design  problems  in  two  domains:   frequency  domain  (ac  analysis) 
and  time  domain  (transient  analysis) .   Variational  calculus 
was  used  to  derive  expressions  that  enabled  the  calculation 
of  all  parameter  gradients  after  one  analysis  of  the  original 
network  and  one  analysis  of  an  associated  network,  known  as 
the  adjoint  network.   Thus,  regardless  of  the  size  of  the 
network  under  consideration,  all  gradients  can  be  determined 
after  two  network  analyses,  an  impressive  advance  beyond  the 
method  of  perturbation.   This  v/ork  was  later  extended  by 
Director  and  Rohrer^"®  to  include  consideration  of  generalized 
network  elements. 

With  an  efficient  gradient  calculation  method  available, 
many  researchers^" ^ ^  developed  automated  design  procedures. 
Existing  techniques,  such  as  network  analysis,  were  easily 
adapted  to  network  design.   Some  work  was  done  on  ac  and  dc 
procedures,  but  more  effort  was  devoted  to  transient  design. 
The  dominance  of  transient  design  was  due  to  its  ability  to 
handle  nonlinear  digital  circuits  which  have  been  used 
extensively  in  recent  high-speed  computers. 

Linear  Amplifier  Design 

An  automated  procedure  for  linear  cimplifier  design  will 
be  described  which  utilizes  the  efficient  "adjoint  netv/ork" 


method  for  gradient  calculation.   This  design  procedure 
differs  from  previous  work  because  it  simultaneously  con- 
siders the  ac  and  dc  characteristics  of  a  circuit. 

A  basic  premise  upon  which  we  base  simultaneous  ac  and 
dc  design  is  that  the  total  operation  of  a  linear  amplifier 
can  be  adequately  represented  by  the  superposition  of  its  ac 
and  dc  responses.   V/hen  only  static  inputs  are  applied  to  a 
linear  amplifier,  dc  analysis  provides  a  complete  description 
of  the  amplifier's  operation.   The  dynamic  behavior  of  a 
linear  amplifier  can  be  ascertained  through  ac  analysis  if 
two  prerequisites  are  satisfied:   the  static  operation  of  the 
amplifier  must  be  known  so  that  ac  models  can  be  specified 
and  all  dynamic  excitations  must  consist  of  sums  of  sinusoids. 
These  dynamic  excitations  must  also  cause  responses  which  are 
small  compared  to  existing  dc  conditions.   Since  satisfaction 
of  these  conditions  is  common,  the  total  operation  of  a  linear 
amplifier  can  be  obtained  through  the  superposition  of  ac  and 
dc  analyses. 

Resistors,  capacitors,  diodes  and  bipolar  transistors 
are  the  most  common  circuit  elements  used  in  integrated  cir- 
cuits.  Resistors  and  capacitors  are  naturally  linear  and 
therefore  may  be  represented  by  linear  models.   But  diodes 
and  transistors,  especially  for  dc  analysis,  are  nonlinear. 
An  iterative  procedure  is  used  to  update  a  linearization  of 
the  networl:  equations  so  that  conventional  linear  solution 
techniques  can  be  employed. 


A  performance  function  is  developed  that  is  capable  of 
representing  all  of  the  ac  and  dc  specifications  with  a 
scalar  measure.   Flexibility  is  built  into  the  function  so 
that  many  different  criteria  may  be  specified.   Typically, 
the  design  specifications  include  requirements  on  amplifier 
gain  and  input  and  output  impedance  over  a  frequency  range 
and  requirements  on  power  dissipation  and  temperature 
sensitivity  over  a  temperature  range. 

Gradient  information  is  determined  by  use  of  the  adjoint 
network  method.   All  parameter  gradients  can  be  obtained  with 
only  two  network  analyses  at  each  frequency  and  two  analyses 
at  each  temperature.   Element  growth  and  removal  have  been 
shown  to  be  feasible  for  resistors  and  capacitors  with  these 
gradients.   Growth  and  removal  of  transistors  are  explored  in 
order  to  allow  the  algorithm  total  flexibility. 

Parameter  adjustment  is  effected  by  a  sophisticated 
algorithm  which  utilizes  the  gradients  to  optimize  the  per- 
formance function.   In  order  to  obtain  realizable  designs, 
the  network  parameters  are  constrained  within  reasonable 
boundaries.   The  presence  of  these  constraints  is  considered 
when  the  optimization  algorithm  is  selected. 

The  main  objective  of  this  research  has  been  to  develop 
an  automated  design  algorithm  that  simultaneously  meets  ac  and 
dc  specifications.  The  computer  program  FROLIC  (Fabrication 
Related  Optimization  of  Linear  Integrated  Circuits)  has  been 
written  and  will  be  used  to  demonstrate  the  effectiveness  of 
the  algorithm.   The  techniques  that  are  incorporated  in  the 
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algorithm  will  be  described  in  this  dissertation.   Chapter  2 
describes  the  models  of  circuit  elements  and  the  method  of 
analysis.   Chapter  3  discusses  the  performance  function,  the 
optimization  strategy,  the  derivation  of  the  gradient  ex- 
pressions, and  the  numerical  integration  procedure.   Chapter 
4  develops  indicators  that  are  used  to  suggest  transistor 
growth  and  removal.   Chapter  5  describes  the  overall  design 
algorithm  and  som>e  methods  that  are  necessary  to  link  together 
the  ac  and  dc  segments.   Chapter  6  contains  some  concluding 
remarks  and  makes  some  suggestions  for  further  research. 


CHAPTER  2 
MODELS  AND  ANALYSIS  TECHNIQUES 

A  basic  part  of  an  automated  network,  design  package  is 
the  analysis  subroutine  which  is  used  repetitively  to  assess 
network  performance.   Since  analysis  is  performed  often,  it 
must  be  accomplished  efficiently.   We  start  with  a  discussion 
of  the  physical  elements  which  can  be  fabricated  in  a  mono- 
lithic integrated  circuit  and  the  ideal  elements  used  to  model 
them.   Tl^en  we  proceed  to  a  discussion  of  equation  formulation 
and  solution  techniques. 

Modeling  of  Physical  Elements 

It  is  important  to  distinguish  between  the  physical 
elements  present  in  an  integrated  circuit  and  the  ideal 
elements  which  are  used  to  model  them.   We  consider  networks 
which  are  made  up  of  integrated  resistors,  capacitors,  diodes, 
and  bipolar  transistors.   Each  of  these  physical  elements 
is  modeled  by  some  of  the  following  ideal  elements:   resistors, 
capacitors,  inductors,  voltage-controlled  current  sources, 
diodes,  and  independent  voltage  and  current  sources.   The 
ideal  element  values  are  specified  by  parameters  which  are 
related  to  the  characteristics  of  the  piiysical  elements,  i.e., 
the  planar  dimensions  of  the  masks,  the  fabrication  procedure, 
and  possibly  the  dc  operating  point.   A  linear  temperature 
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dependence  is  associated  with  each  model  to  allow  for  the 
effect  of  temperature  on  the  performance  of  each  physical 

element. 

The  ac  steady-state,  or  frequency  domain,  behavior  of 
diodes  and  transistors  depends  upon  the  dc ,  or  large  signal, 
operating  conditions  of  the  network.   Therefore,  before  the 
frequency  response  of  an  integrated  circuit  can  be  deter- 
mined, a  dc  analysis  must  be  performed  to  ascertain  the 
operating  points.   The  models  used  to  represent  each  element 
may  be  different  under  dc  and  ac  conditions.   This  fact  is 
discussed  below  in  relation  to  each  element  type. 

Integrated  resistors  are  typically  realized  as  base 
diffusions.   The  v-i  relationship  for  a  resistor  is  approxi- 
m.ated  by  ; 

v^(T)  =  Ps[i^/w]{l+af,[T-T^,]}if,(T) 

where  T  is  temperature  in  degrees  Kelvin, 

V  (T)  is  the  branch  voltage,  -   — 

R 

p   is  the  sheet  resistivity, 
s 

£/w  is  the  length  to  width  ratio  of  the  base  diffusion, 

a   is  a  temperature  coefficient, 
R 

T   is  298"  Kelvin,  and 
0 

i  (T)  is  the  branch  current. 
R 
For  a  given  temperature,  an  integrated  resistor  can  be 

modeled  by  an  ideal  resistor  with  resistance 

R  =  Pg[e/w]{l+a^[T-T^]}. 


Since  this  expression  is  independent  of  frequency  and  oper- 
ating point,  it  is  used  for  both  dc  and  ac  analysis. 

Integrated  capacitors  are  realized  as  reverse-biased 
base-collector  or  base-emitter  junctions  or  metal-oxide- 
silicon  structures  whose  area  A  is  proportional  to  the 
capacitance  of  the  element.   The  i-v  relationship  of  an 
integrated  capacitor  is  approxim.ated  by 

i  (t,T)  =  eA{H-a  [T-T  ]}9v  (t,T)/8t 

where  t  is  time  and  6  is  a  processing  parameter.   For  a 
given  temperature,  an  integrated  capacitor  can  be  modeled  by 
an  ideal  capacitor  with  capacitance 


C  =  eA{l+a  [T-T  1 } 
C      0 


Under  dc  conditions, 


9y  (t,T)/3t  =  0. 

Therefore,  for  dc  analysis 

i^(t,T)  =  0 

and  the  capacitor  may  be  considered  to  be  a  zero-valued 
current  source  or  an  open  circuit.   For  sinusoidal  steady- 
state  analysis,  we  employ  a  complex-valued  branch  relation- 
ship 


I  (o))  --^    juCl/  (to) 
c         c 

where  w  is  the  angular  frequency  of  the  excitation  and  I  (co) 
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and  V    (co)  are  coiriplex-valued  branch  current  and  voltage 
phasors,  respectively. 

Integrated  diodes  are  described  by  the  nonlinear  V-I 
relationship 

V  (T)  =  [nkT/q]£n{l  (T)/[J  A (1+a  (T-T  ) ) ] +1}+I  (T) R 
D  D        S       D      0  D      S 

where  n  is  the  emission  coefficient, 

k  is  Boltzmann's  constant, 

q  is  the  magnitude  of  an  electronic  charge, 

J   is  a  processing  dependent  current  density,  and 

R   is  a  series  resistance. 
S 

The  dynamic  effects  that  influence . the  operation  of  a  diode 
will  be  considered  later.   For  a  given  temperature,  an 
integrated  diode  is  represented  by  an  ideal  diode  whose 
relationship  is 


V   =  [nkT/q]£n[I  /I  +1] , 

D  D   S 


or  its  inverse 


1=1  {cxp[qV  /(nkT)]-l}, 
D      S         D 


in  series  with  resistance  R  .   The  saturation  current  I   is 

s  S 

defined  as  ^       -      .       ■ 


I^  =  J3A[l+a^(T-T^)]. 

During  ac  analysis,  a  diode  is  replaced  by  a  linear  RC  net- 
work whose  dominant  low-frequency  element  i.s  the  small  signal 
resistance  of  the  diode  about  the  dc  operating  point: 
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R  =  9V  /3I  I    =  nkT/(qI  )£Xp[-qV  /  nkT  ] 
D    D  V  S    '      D 

D 

where  V   is  the  dc  junction  voltage.   Two  capacitances  are 
added  to  the  ac  diode  model  to  represent  the  dominant  dy- 
namic effects  in  an  integrated  diode:   a  diffusion  capacitance 
to  represent  the  changes  in  minority  carriers  stored  in  the 
quasi-neutral  regions  and  a  transition  capacitance  to  repre- 
sent the  changes  in  charge  stored  in  the  space-charge  region. 
The  diffusion  capacitance  is  defined  by  the  relationship 

Cp  =  ql2/{nkTa)^)cxp[qV|^/(nkT)] 

where  w   is  an  angular  frequency  processing  constant.   The 
transition  capacitance  is  defined  by  the  relationship 

C^  =  cy(l-Vp/VQ)**(l/m) 

where  C   is  the  transition  capacitance  at  V  =0,  V   is  the 
0  DO 

contact  potential,  and  m  is  the  grading  coefficient.   The 
complete  ac  diode  model  consists  of  the  small  signal  resist- 
ance in  parallel  v;ith  two  capacitances,  all  of  which  are  in 

series  with  the  resistance  R  ,  as  shown  in  Fig.  2-1. 

s 

The  integrated  circuits  under  consideration  employ  both 
NPN  and  PNP  three-layer  bipolar  transistors.   The  following 
discussion  describes  only  the  NPN  transistor  models,  but  the 
PNP  models  are  developed  similarly. 

Separate  consideration  of  the  ac  and  dc  performance  of 
a  transistor  with  distinct  models  for  each  mode  of  operation  . 
allows  simpler  models  to  be  used  without  sacrificing  accuracy. 
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Fig.  2-1 
Model  used  to  represent  a  diode  for  ac  analysis 
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Since  it  is  desirable  to  use  the  same  network,  configuration 
for  each  mode  of  operation,  all  transistor  models  are  re- 
formulated as  two-port  networks  represented  by  y-parameters , 
as  shown  in  Fig.  2-2.   This  formulation  allows  all  transistor 
models  to  be  represented  by  two  branches,  tv70  voltage- 
controlled  current  sources,  and  for  the  dc  case,  two  inde- 
pendent sources.   (Note  that  no  internal  nodes  are  present 
in  this  model. ) 

The  transistor  model  can  be  separated  into  intrinsic 
and  extrinsic  parts.   Extrinsic  effects  are  due  to  the  bulk 
material  between  the  intrinsic  transistor  and  the  emitter, 
base,  and  collector  contacts  and  are  modeled  by  the  three- 
resistors  shown  in  Fig.  2-3.   The  remainder  of  this  discussion 
is  concerned  only  with  the  intrinsic  transistor. 

The  dc  or  large  signal,  performance  of  a  transistor  may 
be  adequately  represented  by  the  Ebers-Moll  model,^'*^^  as 
shown  in  Fig.  2-4,  when  exposed  to  typical  operating  conditions. 
Each  diode  in  this  model  is  assumed  to  obey  the  standard 
diode  relation 


1  =  1. 


X 


{cxpIqVj^^/(n^kT)]-l}  (2.1) 


where  I^ ,  V^^ ,  and  n^  are  I^_'  ^r^i  '  ^^^  ^i  f°^  the  emitter- 

X     DX  X         t  S     L)  1  i 

base  diode  and  I   ,  V^„,  and  n^  for  the  collector-base  diode. 

cs '   D2 '       2 

Temperature  dependence  is  incorporated  in  this  model  through 
the  following  relationships: 


(a) 
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(c) 


Fig.  2-2 

Each  transistor  is  described  as  a  two-port 

(a)  Terminal  voltages  and  currents 

(b)  Y-pararneter  description 

(c)  Circuit  description 
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Ebers-Moll  model  used  to  model 
the  dc  behavior  of  a  NPN  transistor 
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^cs  =  ^csoI^-^^z^T-'^o^J' 


ana 


3p  =  Bpo[l+a3(T-To)] 

where  I   .  r  "'"rqo'  ^""^  ^f   ^""^^  parameter,  values  at  T   and  a  , 
ttg  and  a^  are  temperature  coefficients.   Since  3p  is  not 
directly  used  in  the  Ebers-Moll  model,  the  following  relation- 
ship 


defines  a  which  is  a  controlled-source  gain  coefficient.   It 
F 

should  be  noted  that  dynamic  effects  have  been  omitted  in  the 
above  characterization. 

A  model  that  represents  the  ac,  or  small  "signal,  per- 
formance of  a  transistor  can  be  obtained  directly  from  the 
Ebers-Moll  model  by  linearizing  the  diodes  about  their  dc 
operating  points  and  adding  components  to  represent  the 
dynamic  characteristics  of  each  diode.   Although  this  model 
would  be  suitable,  a  more  convenient  model  can  be  developed 
through  a  simple  transformation.^^   The  first  step  of  the 
transformation  is  accomplished  by  substituting  two  equal 
current  sources  for  each  of  the  two  sources  already  in  the 
Ebers-Moll  model,  as  shown  in  Fig.  2-5 (a) ,  so  that  the  currents 
flowing  in  the  model  are  unchanged.   The  second  step  is  to 
replace  each  parallel  combination  of  a  diode  and  current 
source  with  a  diode  that  has  the  I   of  (2.1)  defined  by 

At 
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Fig.  2-5 

Transformation  of  Ebers-Moll  model 

(a)  First  step 

(b)  Second  step 

(c)  Final  version 


as 


^ES  =  ^ES'^-«F> 


for  the  emitter-base  combination  and 


^^s  =  ^cs(I--r) 


for  the  collector-base  combination.   If  we  define 


l[    =   1^  (l-a^) 


and 


i;  =  i^^i-^R^' 


then  we  have  the  network  of  Fig,  2-5 (b).   By  defining 

Bp  =  Op/d-ap) 


and 


^R  =  «r/(1~«r)' 

we  have  the  final  version  of  the  transformed  Ebers-Moll  model 
shown  in  Fig.  2-5(c). 

To  obtain  a  small  signal  model  from  this  transformed 
model,  each  diode  must  be  replaced  by  a  small  signal  resist- 
ance and  two  capacitors,  as  shown  in  Fig.  2-6.   The  diodes' 
small  signal  resistances  are 

Rj  =  [njkT6p/(apqIp3)]exp[-qV^^/(njkT)J 

and 
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B 


Fig.  2-6 

Linearized  transformed  Ebers-Mcll  model 
used  to  model  the  small-signal  behavior  of  a  transistor 
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Tw^o  capacitors  for  each  diode  represent  the  change  in  minority 
carriers  stored  in  the  quasi-neutral  regions  and  the  change 
in  charge  stored  in  the  space-charge  region.   I'he  diffusion 
capacitances  which  model  the  quasi-neutral  regions  are  defined 
by  the  relations 

and 

where  w  ,  and  w^g  ^^^  angular  frequency  parameters  which  are 
dependent  upon  the  processing  and  the  planar  dimensions  of 
the  masks.   The  transition  capacitances  which  model  the  space- 
charge  layer  are  defined  by  the  relations 

C    =  c   /(1-V   A/   )**(l/m  ) 

and 

C    =  C   /(1-V   /V   )**  (1/m  ) 
T2      02^      D2^  02       '     2 

where  C    and  C    are  the  transition  capacitances  measured  at 
V   =0  and  ^^2!^^,    respectively,  V^  ^  and  V^^  a^^e  contact 
potentials,  and  m   and  m   are  the  grading  coefficients. 

The  linearized  transformed  Ebers-Moll  model  is  a  suit- 
able small-signal  transistor  model  for  use  in  al]  four  regions 
of  transistor  operation.   Two  problems  with  this  model  are 
the  large  number  of  parameters  that  need  to  be  specified  and 
the  absence  of  consideration  of  the  Early  effect.^" 
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For  transistors  knovm  to  operate  in  the  forward- active 
region,  a  simpler  small-signal  model,  the  hybrid-pi  model^^ 
shown  in  Fig.  2-7,  is  available.   This  model  r-equires  fewer 
parameters  and  includes  consideration  of  the  Early  effect. 
This  model  should  not  be  used  for  any  transistor  which  m.ay 
operate  in  a  region  other  than  forward-active,  as  inaccurate 
results  will  be  obtained. 

Nonlinear  DC  Analysis  Procedure 

The  presence  of  diodes  in  a  network  requires  the  solution 
of  nonlinear  equations  in  order  to  obtain  large  signal  re- 
sponses.  Since  efficient  methods  for  the  solution  of  linear 
equations  are  available,  we  are  motivated  to  employ  an  analy- 
sis scheme  which  is  based  upon  repeated  linearizations  of  the 
nonlinear  branch  relationships  about  estimates  of  their 
operating  points.   In  effect,  what  results  is  a  technique 
which  replaces  each  diode  with  a  series  combination  of  a 
conductor  and  an  independent  current  source. 

The  nonlinear  equations  are  of  the  form 

f (V)  =  0  (2.2) 

where  each  equation  is  the  sura  of  currents  leaving  a  node  of 
a  network  and  V  is  the  vector  of  node  voltages.   An  equation 
and  a  voltage  are  associated  v;ith  each  node  of  a  network 
except  one  node  which  is  specified  as  a  reference  or  ground. 
A  suitable  solution  technique  is  the  iterative  Newton- 
Raphson  algorithm  which  can  be  derived  by  expanding  f (V)  in 
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The   hybrid-pi   model 
used    to  model    the    small-signal   beliavior   of    transistors 
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a  Taylor's  series: 

f(V^+l)  =  f  0/^+J(V^)  [V^'+l-V^^]  +  ....  (2.3) 

• 

where  v'"^"^^  is  the  node  voltage  vector  at  iteration  M+i,  V 
is  the  node  voltage  vector  at  iteration  M,  J(V)  is  the 
Jacobian  of  f^(V)  ,  and  higher  order  terms  are  neglected.   If 
V*^"^^  is  assumed  to  satisfy  (2,2),  then  (2.3)  can  be  rev/ritten 
in  the  form 

v^+i  =   v^-[j(v^]-^f(v^.  - -  -  .   - 

This  expression  allows  an  initial  estimate  of  the  node  voltages 
to  be  improved,  iteratively  until  the  actual  node  voltages  are 
obtained. 

An  equivalent  procedure  will  now  be  described  which  allows 
the  nonlinear  equations  to  be  solved  by  repetitive  analysis 
of  a  linearized  network  in  which  each  nonlinearity  is  replaced 
by  its  Taylor's  series  expansion.   The  diode  relationship  is 

and  its  Taylor's  series  expansion  is 

C^  =  i3[c.xp(v'^A^)-l]+i3A^  ^^P(v>^)(v^^^-v^)  +  ... 

where 


V   =  nkT/q 
X        ^ 


If  we  define 


^D  ==  ^s'^V;  cxp(V^J,/V^), 
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then 


I    =  I  [£Xp(V  A  )-i]+G  (v"   -V  )  . 

D        S        D   X        D   D     D 


This  expression  can  be  put  in  the  form 


C'=G»vr'-a3[exp(v>,,-l]-GX) 


which  is  the  I-V  relationship  of  a  current  source  of  value 


I   =1  [£Xp(V  /V  )-l]-G  V 
D      S     '    D^  X     -^    D  D 


in  parallel  with  a  conductance  G  .   This  configuration, 
shown  in  Fig.  2-8,  replaces  each  diode  in  the  network  and, 
upon  repeated  linear  analyses,  yields  the  actual  operating 
points  of  the  network. 

An  example  will  now  show  the  equivalence  of  these  two 
procedures.   For  the  network  in  Fig.  2-9  (a)  ,  we  write  equa- 
tions of  the  form  (2.2)  which  are 


f  (V)  =  0  =  -I+G  V  +1^  JaxpCV  A^)-l] 

1—  AlSl        IX 


-l^^[^Xp({V^-V^)/V^)-l] 


(2.4) 


and 


f^(v)   =  0  =  i^^[e.xp{{v^-v^)/v^)-l]+G^\r^ 


(2.5) 


The  J'acobian  in  (2.3)  is  defined  as 


sfj/avj    ^fj/^v^ 
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Fig.  2-8 

Linearized  diode  model  on  right  replaces  all  diodes 
in  networks  to  be  analyzed 
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(a)  Simple  diode  network  used  in  example 

(b)  Linearized  version  of  the  same  network 


28 


where 


'Dl  =  Isi/^X"P^^>X) 


and 


,M  „M, 


We 


^02  =  i52Axexp((V2-v^)A,,). 
now  put  (2.4)  and  (2.5)  into  the  linearized  form  of  (2.3) 
0  =  -l+G^v'J-M3,[exp(V>,)-l]-Is2[«P((^^^^i^/^^x)"^l 


M    M 


-G 


M 
02 


1      1 


2      2 


and 


o  =  I  Jexp((v^-v';')/v^)-l]+G^v^ 


+ 


-G 


M 
D2 


M 
G   +G 
L  D2   bJ 


^1      1 


M+1   M 
V    -V 
2      2J 


Where  the  superscript  T  indicates  transposition.   The  vectors 
are  multiplied  and  terms  are  grouped  to  form  the  following 
expressions : 


0  =  -^+^Dr'D2 


,(r    +r^    +r^'^    WM  +  I-G^'l  vf^+1 


(2.6) 


and 
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°  -  ^D2  ^D2^1   +^^D2+^B^^2 


(2.7) 


where 


Dl 


M  „M 


,lexp(v';>,)-i)-G'^,v;' 


and 


'^M  MM  M     M   M 

I    =  I    £Xp((V  -V  )A  )-l]-G   (V  -V  ) 

D2      52'^   '2   1  '^  X    -'   D2   2   i 


A  matrix  form  of  (2.6)  and  (2.7)  is 


M    M 
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M    M 
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This  expression  could  have  been  obtained  directly  from  the 
network  in  Fig.  2-9  (b)  v/here  each  diode  is  replaced  by  its 
linearized  model  of  Fig.  2-8.   The  node  voltages  obtained 
from  the  solution  of  these  equations  allov;  determination  of 
new  linearized  models  and  the  process  is  repeated  until  the 
actual  operating  points  are  obtained. 

The  two  Newton-Raphson  procedures  that. have  been  dis- 
cussed will  produce  the  same  result  in  the  same  number  of 
iterations.   The  linearized  model  scheme  is  used  in  this 
algorithm. 

A  modification  to  the  Newton-Raphson  algoritlim  is 
implemented  to  avoid  numerical  problems  and  speed  convergence, 
The  linearized  model  requires  the  evaluation  of  an  exponen- 
tial function  of  the  diode  voltage.   If  the  diode  voltage  is 


large  compared  to  V  ,  the  exponential  function  could  be 
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prohibitively  large  and  cause  numerical  problems  or  possible 
divergence  by  the  iterative  procedure.   Therefore,  at  all 
iterations  where  the  diode  voltage  increases  with  respect  to 
the  previous  iteration,  an  alternate  model  is  used  which  is 
based  on  the  current  that  was  flowing  in  the  linearized  diode 
model.   This  alternate  model  is  defined  by  the  relations 

M       M 
G   =  (I  +1  )/V 
D       D   S^  X 

and 

"M      MM,    „   ,  M  ,    ,  ,  , 
I   =  I  -G  [V  ln{l    /I  +1  ] . 
D      D   d"^  X      D'^  S     •■ 

At  iterations  where  the  diode  voltage  decreases,  the  standard 
Nev/ton-Raphson  method  is  used.   The  two  estimated  operating 
points  about  which  the  diode  model  is  linearized  are  shown  in 
Fig.  2-10.   This  modified  Newton-Raphson  algorithm  has  been 
found  to  converge  reliably. 

Another  numerical  problem  that  is  encountered  during  use 
of  the  linearized  model  is  caused  by  the  wide  range  of  values 
that  can  be  assigned  to  the  resistance.   This  problem  is 
avoided  by  linearizing  the  diode  curve  to  the  right  of  the 
point  v/here  the  resistance  is  equal  to  a  defined  short  circuit 
and  to  the  left  of  the  origin,  as  shown  in  Fig.  2-11. 

The  number  of  linear  analyses  needed  to  obtain  an  accu- 
rate analysis  of  a  nonlinear  network  can  be  reduced  by  the 
use  of  passive  initial  models  for  the  nonlinear  elements.   A 
satisfactory  initial  model  for  non-transistor-model  diodes 
and  emitter-base  diodes  is  a  one-kilohm  resistor.   This 
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Two  possible  estimated  operating  points 
about  which  to  linearize  the  diode  relationship 
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initial  model  assumes  that,  the  diode  is  forward-biased,  but 
if  the  diode  is  reverse-biased,  the  model  will  be  corrected 
after  the  first  linear  analysis.   Collector-base  diodes  are 
assumed  to  be  reverse-biased  and  their  initial  model  is  a 
resistor  with  the  value  defined  above  for  reverse-biased 
diodes.   These  passive  initial  models  are  superior  to  initial 
models  v/ith  current  sources  because  of  the  possibility  of 
initial  divergence  with  the  active  model  when  a  wrong  assump- 
tion is  made. 

Many  criteria  could  be  used  to  determine  when  the  linear- 
ized network  has  accurately  reproduced  the  voltages  and  currents 
of  the  nonlinear  network.   A  possible  choice,  to  specify  a 
maximum  change  allowable  in  node  voltages  at  two  consecutive 
iterations,  would  be  satisfactory  for  some  networks.   Problems 
with  this  criterion  are  the  requirement  that  the  user  know  the 
range  of  approximate  magnitudes  of  the  node  voltages  in 
advance  and  each  node  voltage  will  have  a  different  required 
accuracy  that  is  dependent  on  the  magnitude  of  the  voltage. 
An  improved  method  is  to  determine  the  relative  change  in 
node  voltages  between  iterations  by  the  relation 

I  (v^'^^-V^)/v'-'*'^  I  <  e  (2.8) 

where  vj;:"^"^  is  the  voltage  of  node  N  at  iteration  L+l,  V^  is 
the  voltage  of  node  N  at  iteration  i. ,  and  c  is  a  parameter 
specified  by  the  user.   If  any  node  voltage  fails  to  satisfy 
this  criterion,  another  iteration  should  be  started  unless  a 
maximum  number  of  iterations  has  been  reached.   This  criterion 
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L+  1 

is  sufficient  for  almost  all  netv/orks.   When  V   =0,  an 

N 
alternate  criteria,  |v  |<e,  is  used  which  avoids  division  by 

zero. 

Some  networks  contain  transistors  or  diodes  whose  linear- 
ized models,  during  the  iterative  process,  are  biased  so  that 
their  approach  to  the  true  operating  point  is  extremely  slow. 
Unless  the  value  of  c    is  very  small,  the  node  voltages  change 
so  slowly  that  termination  criterion  (2.8)  could  be  satisfied 
before  the  true  operating  conditions  are  determined.   Such  a 
situation  is  demonstrated  by  the  network  of  Fig.  2-12.   The 
node  voltages  for  each  iteration  of  the  analysis  scheme,  as 
obtained  by  the  computer  program  FROLIC  whose  flow  charts  are 
shown  in  the  appendix,  are  listed  in  Table  2-1.   The  value 
of  e  for  this  analysis  was  l.E-4,  a  value  suitable  for  most 
problems.   The  node  voltages  at  iterations  five  and  six, 
identical  for  five  digits  of  accuracy,  would  have  satisfied 
the  criteria  in  (2.8)  at  iteration  six  and  yielded  an  incorrect 
voltage  for  node  twelve.   This  problem  could  be  alleviated  by 
reducing  the  value  of  c  which  would  have  forced  the  analysis 
beyond  iteration  twelve.   A  more  satisfactory  solution  is  to 
require  all  diode  relationships  to  be  satisfied  within  ten 
percent  before  node  voltages  are  checked,  forcing  the  network 
to  the  neighborhood  of  the  final  solution  before  it  is  poss- 
ible to  terminate  analysis  and  yet  not  overriding  the  flex- 
ibility of  the  check  on  the  node  voltages.   This  combination 
approach  is  used  in  P'ROLIC  and  was  in  effect  during  the  deter- 
mination of  the  node  voltages  in  Table  2-1.   When  the  test 
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Fig.  2-12 
Ten-transistor  differential  amplifier 
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for  satisfaction  of  the  diode  relationship  is  the  only  termi- 
nation criterion,  the  agreement  necessary  in  some  diodes  is 
too  severe  for  others,  causing  prolonged  analysis  in  many 
networks,  while  the  node  voltages  do  not  change  significantly. 

Solution  of  Linear  Network  Equations 

It  has  been  shown  that  nonlinear  dc  network  equations 
can  be  solved  by  repeated  analysis  of  a  sequence  of  linear- 
ized networks.   In  addition,  the  ac  models  which  were  described 
in  a  previous  section  have  linear  1-U    relationships  and  there- 
fore ac  networks  are  linear.   Nodal  analysis  is  employed  to 
generate  linear  equations  which  are  solved  by  sparse  matrix 
techniques. 

The  ac  and  dc  linearized  networks  could  be  analyzed  with 
the  same  set  of  computer  instructions  because  the  general 
procedure  is  identical  for  both  cases.   However,  equations 
describing  the  small  signal  ac  behavior  have  complex  coeffi- 
cients v;hereas  the  equations  describing  dc  behavior  have 
real  coefficients.   Since  complex  operations  consume  more 
time,  separate  solution  siobroutines  are  desirable  for  ac  and 
dc  analysis.   Although  separate  subroutines  are  used,  the 
techniques  employed  by  them  are  the  same  and  will  not  be 
distinguished  in  the  discussion  that  follov/s. 

To  avoid  certain  manipulative  difficulties,  we  assume 
each  branch  to  be  of  the  form  shown  in  Fig.  2-13  and  described 
by  the  relationship 

I  -^    V    iV    -V      )^-I   '-'/   (l/'  -1/   ) 

K      K   K   SK     SK   CK   J   SJ 


38 


n. 
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Fig.  2-13 

General  branch  used  in  the  design  algor j thm 
where  (^j-^^gj)  are  voltages  from  another  branch 
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where  I    is  the  branch  current, 

K 

V  is  the  branch  self-admittance, 
K 

\J        is  the  branch  voltage, 
K 

\J         is  the  source  voltage, 
SK 

I    is  the  source  current, 
SK 

V  is  a  transadmittance , 
CK 

l^    is  the  branch  voltage  of  branch  J,  Jj^K ,  and 

V  is  the  source  voltage  of  branch  j. 
SJ 

All  of  the  individual  branch  relations  may  be  combined  into 
the  matrix  relation 


B      B   B   SB     SB   C   B   SB 

We  may  rewrite  this  expression  as 

I_  =  (Z  +Z  Jii  +1   -(}^+Z)^   ••  (2.9) 

B       B   C   B   SB    B   C   SB 

Kirchhoff  s  current  law  can  be  written  in  the  matrix  fom:i 

AT   =  0  (2.10) 


where  A  is  the  nodal  incidence  matrix.   Multiplication  of 
(2.9)  by  A  and  substitution  of  (2,10)  yields  the  relationship 

A(V;_  +V_   )U   =  A(V  +V_   )U   -Aj_   .  (2.11) 

B   C   B        B   C   SB    SB 

Kirchhoff 's  voltage  law  can  be  written  in  the  matrix  form 

k^\J      =  V  (2.12) 
N    ~B 

where  A^  is  the  transpose  of  the  nodal  incidence  matrix  and 
1/   is  the  vector  of  node  voltages.   Substitution  of  (2.12) 

— N 
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into  (2.11)  yields  the  relationship 


If  we  define 


Aiy_  +V_   )aJv_     =  A(y_  +V_   )U_      -AJ   .  (2.13) 

B   C      N        B   C   SB    SB  • 


y_  =  A{v_  +y_  )aJ 

B   C 


and 


B   C   SB   SB 

then  (2.13)  can  be  written 

'  .  yv      =    T.  (2.14) 

•  •  ^  .         N 

(Notice  that  J_   is  not  used  to  represent  the  identity  matrix 
anywhere  in  Wiis  work . ) 

—   ^Significant  importance  can  be  attributed  to  the  ease  and 
efficiency  with  which  the  nodal  equations  are  generated.   The 
node  numbers  and  general  branch's  component  values  are  suffi- 
cient information  to  form  the  coefficients  that  constitute 
the  nodal  admittance  matrix  y_   and  current  vector  I_.   No  effort 
is  required  to  consider  tree  branches,  links,  or  any  other 
topological  designation,  except  nodes,  in  the  selection  of 
appropriate  unknowns. 

The  accuracy  of  nodal  analysis  is  limited  by  the  accuracy 
of  the  models  used  to  represent  the  physical  elements  and  by 
the  accuracy  of  the  digital  computer  that  is  used  to  solve  the 
nodal  equations.   The  accuracy  of  the  models  can  be  improved 
to  some  extent  if  desired,  but  the  accuracy  of  the  computer 
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-looo" 

1 

"•/' 

= 

1000.01 

1/ 

2_ 

0 

is  fixed  unless  special  programming  techniques  are  employed 
which  may  sever ly  reduce  the  efficiency  of  analysis.   An 
example  of  the  numerical  limitations  of  a  digital  computer 
can  be  shown  for  a  computer  that  is  able  to  retain  three 
digits  of  accuracy  for  each  number  it  stores.   If  the  network 
in  Fig.  2-14  were  to  be  analyzed,  then  the  follov/ing  equations 
must  be  solved: 

"lOOO.Ol 
-1000 

But  the  digital  computer  for  this  example  will  retain  only 
three  digits  of  accuracy,  so  the  equations  will  be  stored  as 

10  00 
-1000 

The  equations  are  no  longer  independent;  the  nodal  admittance 
matrix  is  singular  and  the  equations  can  not  be  solved.   This 
problem  can  be  avoided  by  applying  limits  to  the  network 
parameters  such  as  defined  values  for  a  short  circuit  resist- 
ance and  an  open  circuit  resistance.   The  short  circuit 
resistance  must  be  greater  than  zero  to  avoid  division  by 
zero  when  using  the  general  branch  relationship.   The  range 
of  values,  specified  in  powers  of  ten,  between  a  short  and 
an  open  circuit  must  be  limited  to  the  number  of  digits  of 
accuracy  the  computer  is  capable  of  retaining.   These 
parameter  boundaries  do  not  reduce  the  efficiency  or  appli- 
cability of  nodal  analysis,  but  they  do  eliminate  numerical 
inaccuracy  of  the  type  described  above.   It  is  possible  that 
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Fig.  2-14 
Sample,  network  to  be  analyzed 
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the  boundaries  will  alter  the  value  of  a  parameter,  but  the 
effect  is  usually  negligible. 

A  possible  method  for  the  solution  of  (2. "14)  is  to 
determine  the  inverse  of  the  nodal  admittance  matrix  such 
that 

^     =    V^\.  (2.15) 

N 

The  inversion  of  an  n  by  n  matrix  requires  on  the  order  of  n'^ 
operations  where  an  operation  is  defined  as  a  multiplication 
of  a  division.   Additions  and  subtractions  are  not  counted 
because  the  amount  of  time  required  to  perform  these  opera- 
tions is  insignificant  compared  to  the  time  required  for 
multiplications  and  divisions.   The  matrix-vector  multiplica- 
tion in  (2.15)  requires  n^  operations.   The  total  number  of 
operations  to  solve  (2.14)  by  matrix  inversion  is  approximately 

n^+n^. 

LU    factorization,  a  method  that  produces  a  lower  tri- 
angular matrix,  L_,    and  an  upper  triangular  matrix,  U_,  which 
satisfy  the  relation 

allows  the  solution  of  nodal  equations  in  fewer  operations 
than  matrix  inversion.   This  factorization  is  accomplished 
by  a  Crout  reduction.  ^^   The  L_   matrix  consists  of  zeros  in  the 
upper  right  triangle  and  the  t[  matrix  has  ones  on  the  main 
diagonal  and  zeros  in  the  lower  left  triangle.   If  only  the 
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nonzero  and  nonunity  elements  of  L_   and  U  are  saved,  the  U 
matrix  can  be  superimposed  on  the  L_  matrix  allowing  the  two 
matrices  to  be  stored  in  the  same  area  that  was  occupied  by 
the  original  nodal  admittance  matrix.   Once  the  LU    factor- 
ization has  been  effected,  the  determination  of  V^   is  re- 
duced to  a  forward  substitution  and  a  back  substitution.^ 
The  forward  substitution, 

LX  =  1, 

which  determines  an  intermediate  vector  X  is  followed  by  a 
back  substitution, 

■   ^        ■  til/   =  X, 
N 

to  yield  the  desired  node  voltage  vector.   The  LU    factori- 
zation which  eliminates  the  need  for  the  matrix  inversion 
requires  approximately  n"^/3  operations.   The  forward  and 
back  substitutions  are  performed  in  n   operations.   Conse- 
quently, the  total  operation  count  is  about 

n  /3+n 

and  the  related  execution  time  for  the  solution  of  (2.14)  is 
significantly  reduced  by  the  use  of  LU    factorization,  espe- 
cially for  large  n.   If  additional  sets  of  network  responses 
are  to  be  determined  for  different  excitations,  only  n^  more 
operations  need  be  performed. 

The  nodal  admittance  matrix  of  a  typical  electrical 
network  consists  of  many  zero  coefficients.   If  an  optimum 
order  is  determined  for  the  nodal  equations  that  causes  this 
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sparsity  to  be  maintained  in  the  L_  and  U_   matrices,  many 
operations  can  be  saved  if  operations  involving  multiplication 
by  zero  are  avoided.   An  algorithm' '*' ^  ^  is  used  that  orders 
the  nodal  equations  v/ith  the  only  objective  being  to  retain 
as  much  sparsity  as  possible  in  the  L_  and  t[  matrices.   This 
algorithm  orders  the  nodal  admittance  matrix  coefficients  one 
row  at  a  time  while  checking  ahead  in  a  simulated  LU    factor- 
ization to  determine  which  order  would  introduce  the  fewest 
nev/  nonzero  terms.   A  simple  example  demonstrates  the  advantage 
of  this  ordering  scheme  for  the  network  in  Fig.  2-15  (a).   The 
nodal  admittance  matrix  of  this  network,  shov.m  in  Fig.  2-15  (b)  , 
has  two  zeros  but  its  LU    factorization,  shown  in  Fig.  2-15 (c), 
has  no  zeros  other  than  those  which  are  fixed  by  the  tri- 
angular form  of  the  luatrices  because  of  the  poor  order  of 
the  equations.   The  reordered  nodal  admittance  matrix,  shown 
in  Fig.  2-15  (d)  ,  still  has  the  two  zeros  and  so  does  its  LU 
factorization,  shown  in  Fig.  2-15(e). 

After  the  equations  have  been  ordered,  a  program  SPCRED^^ 
is  used  to  generate  three  subroutines  for  matrices  with 
complex  (ac)  coefficients  and  three  subroutines  for  matrices 
with  real  (dc)  coefficients.   These  subroutines  are  generated 
in  machine  code  that  is  immediately  executable.   Each  set  of 
three  subroutines  has  one  that  does  LU    factorization,  a 
second  for  the  forward  and  back  substitutions  to  determine 
the  solution  to  the  original  nodal  equations,  and  a  third 
does  the  forward  and  back  substitutions  to  determine  the 
solution  to  the  equations  generated  by  transposing  the  nodal 
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l/g      l^  ^^  1^. 
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Fig.  2-15 

(a)  Simple  network 

(b)  Nodal  admittance  matrix  for  network  in  (a) 

(c)  LU  factorization  of  (b) 

(d)  Reordered  nodal  admittance  matrix  of  network  in  (a) 

(e)  LU  factorization  of  (d) 
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admittance  matrix.   This  third  set  of  equations  will  be  used 
for  sensitivity  analysis,  to  be  discussed  in  the  next  chapter. 
The  second  and  third  subroutines  may  be  executed  multiple 
times  for  different  source  vectors . 

The  benefits  of  these  sparse  matrix  programs  are  a 
reduction  in  execution  time  and  roundoff  error  because  of 
fewer  operations  and  reduced  storage  space  because  of  fewer 
nonzero  coefficients.   Execution  tim.e  is  also  decreased 
because  general  expressions  with  variable  subscripts  are  not 
necessary.   The  operands  associated  with  each  operation  are 
specified  by  exact  addresses.   Therefore,  the  time  used  to 
determine  addresses  of  operands  through  arithmetic  operations 
on  subscripts  is  saved.   Retaining  only  nonzero  coefficients 
saves  a  significant  amount  of  storage  for  large  values  of  n. 
Another  feature  of  the  sparse  matrix  technique,  the  variable 
length  of  the  generated  subroutines,  is  dependent  on  the 
sparsity  of  the  nodal  admittance  matrix.   For  a  sparse  matrix, 
the  subroutines  will  have  a  length  comparable  to  that  of  a 
subroutine  written  to  solve  general  sets  of  equations ;  if 
the  matrix  is  not  sparse,  the  three  subroutines  could  get 
prohibitively  large.   Few  networks  are  described  by  equations 
with  a  dense  nodal  admittance  matrix  so  this  problem  is 
seldom  significant. 

Since  the  analysis  subroutine  is  to  be.  used  frequently 
during  design,  we  require  that  the  techniques  which  are 
employed  be  efficient,  accurate,  and  suited  to  the  needs  of 
the  design  algorithm.   Nodal  analysis,  as  will  be  seen  later, 
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meets  these  requirements  v/hen  used  in  conjunction  with  the 
adjoint  network^  method  for  sensitivity  calculation. 


CHAPTER  3 
PERFORMANCE  FUNCTIONS  AND  MINIMIZATION  METHODS 

Automated  network  design  is  based  upon  the  minimization 
of  a  performance  function  that  suitably  reflects  the  design 
specifications.   In  this  chapter  we  discuss  the  least  pth 
performance  function,  the  conjugate  gradients  minimization 
scheme,  and  the  efficient  and  accurate  determination  of  the 
parameter  gradients  by  the  adjoint  network  method. 

Evaluation  of  Network  Performance 

A  perform.ance  function  assigns  a  scalar  measure  to  the 
ability  of  a  network  to  meet  design  specifications.   In 
particular,  the  performance  function  represents  a  comparison 
of  a  network's  actual  responses  with  the  desired  responses. 
For  the  case  of  linear  integrated  circuit  design,  this 
functional  relationship  must  account  for  both  ac  and  dc 
performance.   Flexibility  can  be  built  into  the  function 
through  weights  which  can  emphasize  the  more  important 
parts  of  the  responses. 

Since  ac  and  dc  responses  are  determined  separately, 
the  ac  and  dc  portions  of  the  performance  function  are 
determined  in  conjunction  with  the  appropriate  analyses  and 
the  resulting  error  values  are  sunmed: 
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E  =  E   +E   .  (3.1) 

AC   DC 

Because  of  the  form  of  the  performance  function,  a  compromise 
can  be  reached  between  the  ac  and  dc  design  specifications. 
A  desirable  performance  function  is  based  on  the  difference, 
or  response  error,  between  the  actual  and  the  desired  re- 
sponses.  All  components  of  the  performance  function  should 
be  nonnegative  so  that  some  components  do  not  cancel  others. 
Raising  each  response  error  to  an  even  integer  power  guar- 
antees that  no  negative  components  exist.   A  suitable  ac 

perfoormance  function  is 

n     w 
AC    F  ^  ■ 

E        =      I      j         l/p[W  (oo)  |X  (aj)/X  (a))-l|]^  dco       (3.2) 

^^    J=l  ui  ^  -^  ^ 


where  n    is  the  number  of  ac  ports, 

M  v.. 

CO   is  the  initial  angular  frequency, 
"'   0)   is  the  final  angular  frequency, 
p  is  a  positive  even  integer, 
W   is  the  nonnegative  weight  at  port  J, 
X  ((jj)  is  the  actual  voltage  or  current  response,  and 

/v  * 

X  (w)  is  the  desired  voltage  or  current  response. 

The  magnitude  of  the  response  error  is  raised  to  the  pth  power 

to  avoid  a  complex  performance  function  as  the  responses  are 

complex  phasors .   A  suitable  dc  performance  function  is 

^AC  +  ^DC  F                     /v 

Epc  "■      I  I         l/p{W  (T)  [X  (T)/X  (T)-1]}P  dT    (3.3) 

J=n   +1  T        J      J      J 

AC  I 

where  n    is  the  number  of  dc  ports, 
T   .is  the  initial  temperature. 
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T   is  the  final  temperature, 

X  (T)  is  the  actual  voltage  or  current  response,  and 
X  (T)  is  the  desired  voltage  or  current  -response. 
The  dependence  of  the  performance  function  on  the  netv/ork 
parameters  is  implicit  in  that  the  actual  responses  are 
dependent  on  the  network  parameters.   To  shov;  this  dependence, 
(3.1)  can  be  v/ritten  as 


E(r,to,T)  =  E^^(r,co)+Ej^^(r,T) 


where  r  is  the  vector  of  network  parameters .   The  dependence 
of  the  ac  performance  function  on  temperature  is  ignored 
because  it  is  assumed  that  the  dc  specifications  will  include 
temperature  desensitization  requirements,  thus  minimizing  the 
effect  of  temperature  on  the  network. 

The  flexibility  of  this  performance  function  is  not 
limited  to  the  use  of  weights.   The  power,  p,  to  which  the 
normalized  response  difference  is  raised,  may  be  assigned 
any  positive  even  integer  value.   For  p=2,  the  function 
represents  a  least-squares  criteria.   As  p  is  increased,  the 
larger  response  errors  experience  increased  emphasis,  and 
as  p  becomes  very  large,  the  performance  function  approaches 
a  minimax  criterion^ ^  of  the  form 

E  =  max[W  JX  /X  -l]]  . 
J    J   J   J 

The  normalization  included  in  the  performance  function 
equalizes  the  effect  of  each  response  on  the  total  error. 
For  example,  the  contribution  from  a  response  of  one  milliamp 
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is  comparable  to  the  contribution  from  a  response  of  one  volt. 
Therefore,  responses  with  smaller  magnitudes  will  not  be 
neglected.   To  avoid  division  by  zero,  desired  responses  of 
zero  magnitude  are  not  normalized  and  the  quantity  that  is 
raised  to  the  pth  power  is  of  the  form 

WJXJ.  (3.4) 

This  alternate  form  has  been  found  to  be  adequate  for  the 
few  occasions  where  it  was  necessary. 

Frequently,  only  the  magnitude  of  an  ac  response  is 
specified  and  the  phase  is  ignored.   For  these  cases,  an 
alternate  form  of  the  ac  portion  of  the  performance  function 
is  used: 

E^^  =  I   /    l/p{W  (co)  [[X  fw)  1/|X  (co)  I-1.]}P  dw.  ■    (3.5) 
-   -  '^^   J=l   u       ^      ^       ^ 

As  before,  if  a  magnitude  of  zero  is  desired,  the  alternate 
form  of  (3.4)  is  used. 

The  frequency  response  of  a  network  is  often  specified 
over  several  decades  of  frequency.   So  that' errors  at  the 
lower  frequencies  are  weighted  the  same  as  those  at  the  higher 
frequencies,  a  logarithmic  scaling  of  frequency  can  be 
employed.   The  performance  functions  (3.2)  and  (3.5)  become 

^AC   ^F   . 

E    =^   J    l/p{W  (u)  iX  (u)/X  (u)-lI}P  du 

'^^   J=l   u        ^      ^      ^ 
I 

and 


E^^  =  I      j         l/p{Wj(u)[|Xj(u) |/]X^(u) 1-1]}^  du,      (3.6) 
j=i   u 


respectively,  where 

u  =  log    CO. 

The  use  of  large  values  of  p,  p>6/has  caused  numerical 
overflows  and  underflows  for  problems  that  have  responses 
that,  are  not  close  to  the  desired  responses  at  the  deter- 
mination of  the  initial  error.   This  problem  can  be  avoided 
by  using  a  small  value  of  p  to  obtain  a  partial  design  and 
then  performing  the  remainder  of  the  design  with  a  larger  p. 

Minimization  of  the  Performance  Function 

Systematic  adjustment  of  the  network  pararaeters  by  a 
sophisticated  multidimensional  search  algorithm  minimizes  the 
performance  function  and  yields  a  network  whose  responses 
approximate  the  desired  responses.   A  desirable  feature  of 
such  a  minimization  algorithm  is  that  it  require  as  few 
function  evaluations  as  possible. 

Many  algorithms  are  available  that  effectively  minimize 
general  functions.   Such  algorithms  are  separable  into  two 
types:   gradient  methods,  i.e.,  those  that  use  partial 
derivations  to  locate  a  minimum,  ^  ®~  ^ "  and  direct  methods  .  ^  ^~'' ^ 
A  comparison  of  the  number  of  function  evaluations  needed  to 
minimize  some  test  functions  gives  an  indication  as  to 
whether  or  not  it  is  desirable  to  use  a  minimization  algorithm 
that  requires  partial  derivatives.   A  comparison  has  been 
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made^"  that  shows  the  two  types  of  methods  are  comparable 
if  the  determination  of  each  parameter  gradient  requires  a 
function  evaluation.   Therefore,  n+1  function  evaluations 
are  assumed  necessary  to  evaluate  a  function  of  n  variables 
and  determine  the  associated  first  partial  derivatives.   If 
a  method  was  available  that  required  fewer  function  evalu- 
ations to  obtain  the  parameter  gradients,  the  gradient 
methods  would  be  superior. 

An  efficient  method  of  determining  partial  derivatives 
has  been  developed^  and  through  use  of  this  method,  all  first 
partial  derivatives  of  the  performance  function  with  respect 
to  the  network  parameters  can  be  determined  with  one  function 
evaluation  and  the  simultaneous  analysis  of  a  related  network. 
This  method  will  be  discussed  in  the  next  section. 

Because  partial  derivatives  can  be  computed  so  efficiently, 
gradient  type  minimization  methods  are  employed.   In  particular, 
two  such  methods  are  implemented:   conjugate  gradients    and 
steepest  descent. ^^   Other  methods,  such  as  the  methods  of 
Fletcher  and  Powell, ^^  Jacobson  and  Oksman,^^  and  Fletcher,'^ 
are  available  and  have  been  found  to  be  more  efficient  than 
either  conjugate  gradients  or  steepest  descent  when  used  to 
locate  an  unconstrained  minimum.   These  superior  methods  are 
not  used  here  because  of  the  presence  of  constraints  on  the 
network  parameters.   Although  all  of  the  methods  which  have 
been  mentioned  are  designed  for  minimization  of  functions  of 
unconstrained  parameters,  conjugate  gradients  and  steepest 
descent  methods  adjust  each  parameter  along  a  direction  that 
is  mostly  independent  of.  the  adjustment  of  the  other  parameters. 
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This  independence  prevents  an  accumulation  of  erroneous  data 
which  is  due  to  parameter  changes  that  were  suggested  at 
previous  iterations  and  had  been  ignored  in  oi'der  to  keep 
parameters  within  boundaries.   For  instance,  the  method  of 
Fletcher  and  Powell  updates  the  direction  vector  with  a 
vector  that  is  dependent  upon  all  parameter  gradients  and 
all  parameter  changes.   Once  a  parameter  encounters  a  bound- 
ary, all  components  of  the  direction  vector  become  distorted 
by  the  inconsistancy  caused  by  the  constraint.   The  method 
of  steepest  descent  adjusts  each  parameter  in  the  negative 
gradient  direction  at  all  times  and  thus  allows  no  inter- 
action between  parameters.   The  method  of  conjugate  gradients 
adjusts  parameters  along  a  direction  which  is  the  negative 
of  the  sum  of  the  present  gradient  component  and  a  portion 
of  recent  gradient  components ;   the  interaction  betv/een 
parameters  is  limited  to  the  determination  of  a  number  that 
is  used  to  scale  the  past  gradient  information  before  it  is 
combined  with  the  current  gradient.   The  scale  factor  is 
dependent  on  the  magnitudes  of  the  gradient^  and  not  on  the 
actual  parameter  changes,  thus  avoiding  introduction  of 
significant  errors  when  a  parameter  encounters  a  constraint. 

The  network  parameter  values  must  be  forced  to  remain 
within  specified  boundaries  in  order  to  consistantly  produce 
netv.'ork  designs  that  can  actually  be  realized  with  a  pre- 
defined fabrication  procedure.   In  particular,  if  the  beta 
of  ci  transistor  was  an  unconstrained  parameter  and  its  value 
was  specified  to  be  a  million,  the  design  would  have  to  be 
discarded  because  the  large  beta  was  not  obtainable. 
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The  method  of  steepest  descent  is  based  on  the  adjust- 
ment of  parameters  in  the  direction  of  the  negative  gradient; 


r   =  r-a2_ 


where  r   is  the  vector  of  new  parameter  values, 
r  is  the  vector  of  old  parameter  values , 
a  is  a  variable  of  the  minimization  algorithm, 
2_  is  the  parameter  gradient  vector. 
This  method  is  very  effective  when  the  initial  parameter 
values  are  poorly  chosen,  but  is  often  ineffective  when 
searching  through  a  narrow  valley  of  the  performance  function, 
Information  concerning  previous  parameter  adjustments  or 
previous  gradients  is  not  retained  and  the  method  has  a  ten- 
dency to  oscillate  instead  of  following  the  contours  of  the 
performance  function. 

The  method  of  conjugate  gradients  is  different  from  the 
method  of  steepest  descent  in  that  the  direction  of  search 
is  a  combination  of  past  and  present  gradient  information. 
The  direction  vector  along  v/hich  the  parameters  are  varied 
is  initially  the  negative  of  the  gradient  vector  and  is 
modified  at  each  following  iteration  by  the  formula 


d    =  -2.   +3  d 
l+l      l+l   II 


where 


I    i+i  i+i    II 
d    is  the  direction  vector  at  iteration  i+i  and  d   is  the 
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direction  vector  at  iteration  i .  The  vectors  d  ,  d  ,  . . . ,  d 
are  A-conjugate  where  the  performance  function  is  assumed  to 
be  of  the  form 

E(r)  =  E(r)+i2(r-r)^A(r-r)+  ... 

where  r  is  the  parameter  vector  at  the  minimum  of  the  perfor- 
mance function  and  A  is  a  symmetric  positive  definite  matrix. 
The  minimum  of  such  a  quadratic  function  of  n  variables  will 
be  obtained  in  n+1  iterations  if  no  constraints  are  present. 
Typically,  the  performance  function  we  use  is  not  quadratic 
and  requires  more  than  n+1  iterations.   It  has  been  found 
that  after  n+1  iterations,  if  the  direction  vector  is  again 
set  equal  to  the  negative  of  the  gradient,  faster  convergence 
is  achieved.   There  are  three  major  reasons  for  this  pro- 
cedure:  the  function  is  generally  not  quadratic  and  thus  the 
direction  vectors  are  not  A-conjugate,  erroneous  data  accum- 
ulate due  to  constraints  and  the  accumulation  of  rounding 
errors.   The  direction  vector  is  also  reset  to  the  negative 
gradient  if  the  value  of  the  performance  function  is  not 
reduced  a  specified  percentage. 

The  minimization  strategy  used  in  this  design  algorithm 
is  based  on  the  conversion  of  a  multidimensional  search  into 
a  one-dimensional  search.   The  parameters  are  adjusted  along 
a  direction  vector,  obtained  through  one  of  the  above  methods, 
until  the  value  of  the  performance  function  no  longer  decreases, 
thus  bounding  the  location  of  the  minimum  along  this  direction. 
Cubic  interpolation  is  then  used  to  obtain  the  minimum  from 
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whence  a  new  direction  is  generated.   This  process  is  re- 
peated until  one  of  four  criteria  is  satisfied:   the  value 
of  the  performance  function  is  less  than  desir.ed,  the  value 
of  the  performance  function  is  changing  less  than  desired, 
an  excessive  number  of  iterations  have  been  executed,  or  an 
excessive  amount  of  time  has  been  consumed. 

A  change  of  variable  is  made  to  avoid  the  elongated 
contours  of  the  performance  function  that  arise  when  the 
designable  parameters  have  values  that  differ  by  many  orders 
of  magnitude.   In  particular  we  use 

-  Y   =  £n (r  ),  1=1,2,3, ... ,n  ■ 
I        I 

where  Y   is  the  new  parameter  and  r   is  the  old  parameter. 
The  new  gradients  can  be  determined  directly  from  the  old 
gradients  by  the  relationship 

8E/3Y   =  3E/8r  [8Y  /3r  ]"^  =  [1/r  ]"^3E/3r,  =  r  3E/9r  . 
/j      'i^ll-^  I         I     II 

This  transformation  has  improved  convergence  characteristics 
and  allows  all  parameters  to  be  varied  similiarly  during  the 
minimization  of  the  performance  function. 

Efficient  Determination  of  Parameter  Gradients 

As  indicated  earlier,  we  require  knowledge  of  the  first- 
order  partial  derivatives  of  the  performance  function  with 
respect  to  the  network  parameters.   A  method  proposed  by 
Director  and  Rohrer^  for  evaluation  of  partial  derivatives 
through  the  use  of  an  auxiliary  network  related  to  the 
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original  network,  has  been  shown  to  be  very  efficient  and 
accurate  when  used  for  network  optimization. 

All  the  elements  (resistors,  capacitors,  diodes,  and 
transistors)  of  a  monolithic  integrated  circuit  are  ultimately 
specified  in  terms  of  planar  dimensions  (the  masks)  and  a 
series  of  diffusions,  oxidations,  etc.  (the  processing  steps). 
Typically  for  economic  reasons  the  processing,  or  impurity 
profiles,  are  fixed  and  the  mask  geometries  varied  to  realize 
the  designable  elements.   The  only  constraint  on  the  geometries 
is  dictated  by  the  resolution  capability  of  the  processing. 
In  particular,  a  resistor  with  resistance  R  is  realized  as  a 
base  diffusion  whose  length  to  width  ratio  is  such  that 

R  =  {l/vj)p    [l+a  (T-T  )] 

S      R  .    0 

where  the  sheet  resistivity  and  the  temperature  coefficient 
are  determined  by  the  processing.   Thus  we  can  consider  R  to 
be  the  designable  electrical  parameter  of  a  resistor  and 
length  to  width  ratio  the  corresponding  designable  mask 
parameter.   A  capacitor  is  realized  as  a  reverse  biased  base- 
collector  or  emitter-base  junction  or  a  metal-oxide-silicon 
structure  whose  area  is  such  that 

C  =  0A[l+a  (T-T  ) ] 
C     0 

where  0  and  the  temperature  coefficient  are  processing  de- 
pendent constants.   The  designable  electrical  parameter  of 
a  capacitor  is  C  while  the  corresponding  designable  mask 
parameter  is  area.   A  junction  diode  has  area  such  that 
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I   =  J  A[l+a  (T-T^) ] 
S     S      D     0 

where  the  current  density  and  the  temperature  coefficient  are 
determined  by  the  processing.   Thus  we  can  consider  I^  the 
designable  electrical  parameter  of  a  diode  and  area  the 
corresponding  designable  mask  parameter.   DC  bias  parameters 
will  be  of  importance  in  the  sequel;  for  the  diode,  such  a 
parameter  is  the  junction  voltage. 

It  is  not  so  apparent  what  the  designable  electrical 
parameters  for  a  transistor  are,  but  a  suitable  choice  for 
designable  mask  parameters  is  the  emitter  area  A   and  perim- 
eter P  .   Actually  any  two  independent  specifications  of 
length  and  v/idth  of  the  emitter  suffice;  these  two  prove 
most  convenient.   Once  the  emitter  geometry  is  given,  the 
base  and  collector  geometries  are  determined  by  the  masking 
tolerances.   Both  3   and  oo   are  completely  specified  for  a 
given  A  ,  P  ,  I  ,  and  the  processing  parameters  a   and  b  : 

EEC  I        J 

1/3   =  a  +a  (P  /A  )+a^(I^/A^)+a  (l/I^) 
'  F      1   2   E   E    3   C   E    4      C 


1/03^  =  b^+b^(A^/I^)+b3A^. 

It  is  apparent  that  the  designable  electrical  parameters  for 

a  transistor  are  1/6  ,  l/oo  ,  and  at  least  one  bias  parameter 

to  specify  I  .   For  the  linearized  Ebers-Moll  ac  model,  we 

choose  the  two  diode  voltages  at  T  ,  V    and  V   .   For  the 

0    Dl       D2 

hybrid-pi  ac  model,  we  choose  the  magnitude  of  the  collector 
current  at  T  ,  |l  ] ,  which  is  a  component  of  g  . 
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An  auxiliary  network  will  be  defined  that  is  topolog- 
ically  identical  to  the  original  network.   It  will  be  shown 
that  this  auxiliary  network,  called  the  adjoint  network, 
consists  of  elements  that  are  related  to  the  elements  of  the 
original  network  and  is  excited  by  error  signals  determined 
from  responses  of  the  original  network.   The  nodal  adm.ittance 
matrix  of  the  adjoint  network  has  been  shown  to  be  the 
transpose  of  the  nodal  admittance  matrix  of  the  original 
network. ^^   Since  the  nodal  admittance  matrix  is  generated 
for  the  analysis  of  the  original  network,  the  analysis  of 
the  adjoint  network  reduces  to  a  for>/ard  and.  back  substitution 
using  the  subroutine  described  in  the  previous  chapter  to 
solve  the  equations  for  the  transposed  nodal  admittance  matrix, 
This  feature  of  the  adjoint  network  in  combination  vri  th  the 
use  of  nodal  analysis  makes  this  method  of  gradient  evaluation 
very  efficient. 

Tellegen's  theorem^®  which  is  used  to  develop  the  re- 
lations needed  to  calculate  the  first-order  partial  deriv- 
atives can  be  written 

n^ 

I' 
I  I 


'B 
I      VJ^    =  0  (3.7) 

1  =  1 


and 

"b 

y  J    m      =0  (3.8) 

1=1    '  I 

where  n   is  the  number  of  branches  in  the  original  network, 
B 

[/      is  the  voltage  of  branch  I  in  the  original  network, 
<J)   is  the  current  of  branch  i  in  the  adjoint  network. 
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I   is  the  current  of  branch  I  in  the  original  network, 
I 

1'   is  the  voltage  of  branch  I  in  the  adjoint  network. 
If  the  parameters  of  the  original  network  are  'perturbed, 
Tellegen's  theorem  is  still  valid: 

I       (l/j+Al/j)$j    =    0    .  (3.9) 

1  =  1 

and 

n 
B 

1  =  1 

If  the  difference  between  (3.7)  and  (3.8)  is  subtracted  from 
the  difference  between  (3.9)  and  (3.10),  the  following 
relation  is  obtained: 

II     II 


I       (T  +AT  )^>      =   0.  (3.10) 


B 

I       (Al/^$^-AT^T^)  =  0.  (3.11) 

1  =  1 


Substitution  of  the  branch  relationships  for  each  element 
type  with  perturbed  parameters  into  (3.11)  will  suggest 
gradient  and  adjoint  element  relationships. 

The  branch  relations  that  are  considered  in  the  ac  and 
dc  cases  for  a  resistor  are 


R       R 


and 


where  G=l/R.  The  purpose  of  the  two  branch  representations 
for  a  resistor  will  become  clear  in  the  next  chapter.  Per- 
turbing the  parameters  and  substituting  the  perturbed  branch 
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relations  into  (3.11)  yields' 


^^^-^^R^R  =  (^^R-^r)^^R-^^R^R^^  (3.12) 

for  resistances  and 

Al/^O^-AT  H'.  =  ($^-G^V)Al'^-l/^^VAG.  (3.13) 

for  conductances.   Two  branch  relations  are  also  used  for 
both  capacitances  and  inductances  in  the  ac  case  to  yield 

Al/^$^-AI^H'^  =  ($^-ja)C'i'^)Al'^-jwl/^H'^AC         (3.14) 

for  capacitances , 

Al/g^g-AI^Tg  =  [l/(jco)S<Dg-'F^]A72+l/(joj)l2$^AS     (3.15) 

for  elastances  where  S=l/C, 

Al/  ^    -AI  '1'   =  (juL'J'  -T  )AI  +jwT  $,  AL         (3.16) 

for  conductances,  and 

Al/j,$^-AI^^'^  =  [$p-l/(jao)r>l'^]Al/^-l/(joj)t/^H'^Ar     (3.17) 

for  reciprocal  inductances  where  r=l/L.   In  each  of  the 
above  relations,  dependence  on  hV    and  AI  can  be  eliminated 
by  choosing  the  adjoint  network  element  to  be  identical  to 
the  element  in  the  original  network. 

A  voltage-controlled  current  source  is  represented  by 
the  two-branch  model  shown  in  Fig.  3-1  and  the  branch  relation- 
ships 

T     =  g  1/ 

VDI      M  VCI 
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vci 


VDI 


Fig.  3-1 

Two-branch  .,odel  for  a  voltage-controlled  current 
where  the  subscripts  vci  and  VDI  indicate 

in\he"or?"'"?  """^  -'■'P'^n^'ent  branches 
in  the  originf-.l  network,  respectively 


source 
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and 

« 

for  both  ac  and  dc  cases.   Substituting  these  expressions 
into  (3.11)  yields 


Al/    <!)    -AT    ¥    +AI/    $    -A  I    'i' 
VCI  VCI     VCI  VCI     VDI  VDI     VDI  VDI 


-     (^         -q  Y    )  M/         +<i>         Al/    -V         ^         Aq  ,         (3.18) 
^  VCI  ^M  VDl'   VCI   VDI   VDI   VCI  VDI  ^M  V^.-i-o; 

The  adjoint  netv.'ork  element  is  chosen  to  be  a  voltage-controlled 
current  source  with  the  controlling  and  dependent  branches  re- 
versing rolls  as  indicated  in  the  branch  relations 


(J)     =  g  H' 
VCI      M  VDI 


and 


$    =  0. 

VDI 


The  dc  diode  relation  with  the  value  of  the  saturation 
current  perturbed  can  be  written 

« 

V  +AV   =  V  ln[{I   +AI  )/{!  +AI  )+l]+(I  +AI  )R  . 
D    D     XD    DS    S         D    DS 

Substituting  into  this  expression  the  first-order  terms  of  a 
Taylor's  series  for  the  logarithm  yields 


V^^D  =  ^X^^^^^d/^S+I^+^^d/^V^s) 


-I  AI  /[I  (I  +1  )]}+(I  +AI  )R  , 
D   S^  "^  S   D   S  -^      D    D   S 
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Therefore, 

Substituting  this  expression  into  (3.11)  yields 

which  can  be  rearranged  into  the  form 

{[V  /(I  +1  )+R  ]$  -"¥    }AI  -V  I  <J)  /[I  (I  +1  )]AI  .       (3.19) 
•^  X    D   S    S   D   D    D   X  D  D    S   D   S  "^   S 

At  this  point,  the  goal  is  to  eliminate  dependence  of  this 
expression  on  AI   by  choosing  the  adjoint  branch  relationship 
to  be  .    _ 

~  "  """      y   =  [V  /(I  +1  )+R  1$  . 

D       X    D   S    S   D 

Thus  the  adjoint  netv/ork  element  of  the  dc  diode  relationship 
is  a  linear  resistor  whose  value  is  the  resistance  of  the 
linearized  diode  model  at  its  operating  point  plus  the  series 
resistance.   The  remainder  of  the  expression  will  be  used  to 
determine  the  gradient  later  in  the  section. 

For  ac  analysis  a  diode  is  represented  by  the  combination 
of  two  resistors  and  two  capacitors  shown  in  Fig.  3-2.   We 
will  omit  consideration  of  the  series  resistance  for  this 
derivation  as  it  can  be  represented  as  an  identical  resistcince 
in  the  adjoint  network  and  it  is  independent  of  the  saturation 
current  and  the  junction  voltage  which  are  the  parameters  of 
interest  in  the  diode.   The  branch  relationship  of  the 
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+ 


+ 


Fig.  3-2 
Model  used  to  represent  a  diode  for  ac  analysis 
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remainder  of  the  ac  diode  model  with  the  value  of  the  satu- 
ration current  and  the  junction  voltage  perturbed  is 

^D+^^D  =  ni3+Al3)/V^£Xp[(V^+AV^)/V^][l+ja3/co^] 

+  jwC  [1-(V  +AV  ) /V       ]**  (1/m)  }  (l/'+Al/') 
0       D    D    DO        ■     D    D 

f 

where  \J      is  the  ac  diode  voltage  less  the  ac  voltage  dropped 

across  the  series  resistance.   If  v/e  substitute  the  linear 
terras  of  a  Taylor's  series  where  necessary,  subtract  the 
unperturbed  expression  from  the  perturbed  form,  and  neglect 
second-order  terms,  we  obtain  the  following  expression: 


-jwC  /[mV   (1--V  /V   )**(l+l/ra)j  I/'av 

■^   0  ■•   DO      d'     do         ^    -•   D   D 


+  lAj^c.xp(v^Ax)  [l+ja3/«^]f|^AIg 


+  {i^A^exp(v^A^)  [l+jw/oj^] 


+ja)CQ/(l-V^ADo)**(Vm)}Al/j^ 


Substitution  of  this  expression  into  (3.11)  yields  the 
following: 
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^Vd-^^d^d  =  -^Is/^x^'^p(^d/^x^[^+^'^/'^x^ 


I     1 

•jtaC    /[mV       (1-V    /V      )  **  (1+1/rc)  1  }  i'    H'    AV 
-^        0  DO  D        DO  -^        D    D       D 


-l/v^exp(v^/v^^)  [1+j^>x^^dV^s 


■^^V^s/^X^'^P^^D/^X^tl-^^a./c^]^'^ 


-jwCQ/(l-Vj^ADo)**(Vm)H'j!^}Al/j!,  (3.20) 

where  H'   is  the  ac  diode  voltage  in  the  adjoint  network  less 
the  voltage  dropped  across  the  series  resistance.   The  adjoint 
network  ac  branch  relationship  of  a  diode  is  therefore 

$  =  {I  /v  exp  (V  /V  )  [1-i-jw/w  ]  • 

D       S   X   '    D   X      -^    X 

+  jajC  /(1-V  /V   )**(l/m)}T' 
■^   0       D   DO  D 

which  is  identical  to  the  ac  diode  branch  relationship  in 
the  original  network.   (Remember  that  the  series  resistance 
was  set  aside  as  an  individual  element  and  is  still  present.) 

For  the  determination  of  the  adjoint  network  elements 
and  gradient  expressions  of  the  transistor  models,  vie   will 
consider  only  the  intrinsic  portions  of  the  models.   The 
three  bulk  resistors  that  have  been  designated  the  extrinsic 
elements  have  corresponding  resistors  as  adjoint  elements. 
Since  these  resistors  are  not  dependent  upon  the  parameters 
of  the  transistor,  their  omission  from  gradient  consideration 
is  justified. 
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First  we  consider  the  dc  Ebers-Moll  model  for  an  NPN 
transistor  which  can  be  expressed  by  the  following  equations: 


where 


V   =  n  kT/q 
XI     1   ^  -^ 


and 


V    =  n  kT/q. 
X2      Z       '  ^ 

For  this  model,  we  consider  the  effect  of  the  perturbation 
of  only  1/B   which  can  be  determined  by  perturbing  the  for- 
ward alpha  and  transf orr.ing  the  result.   The  perturbation  of 
the  forward  alpha  and  use  of  the  first-order  terms  of  the 
Taylor's  series  for  the  exponential  yields 


■E--E  -   -Est^^P^-^EB/^Xl^-l-^^P^-^EB/^Xl^^^EB/^Xl] 


I^+AI^  =  -I, 


+a^l^^[exp(~V^3/V^^)-l-exp(-V^^A,2)AV^^A,2]   (3.23) 


and 


I  +AI   -  (a  +Aa  )I  \zxx>{-^     A   )-l-exp(-V  A   )  AV   A   ] 

C    C       F    F   ES"-   '     EB   XI        '     EB^  XI     Eb'  X 1  ' 


-I   [£xp(-V  A   )-l-exp(-V   A   )AV  A   ]•     (3.24) 

CS         CB   X2  CB   X2     CB   X2 

If  we  subtract  (3.21)  from  (3.23),  s\abtract  (3.22)  from  (3.24), 
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and  neglect  second-order  terms,  we  arrive  at  the  following 
expressions: 

AI  =  I   /V   [exp(-V  /V  )]AV  -a  (I  A   )[exp(-V  /V  )]AV 

E     ES   XI    '     EB^  XI  -"   EB   R   CS^  X2  "^   '     CB'  X2  ■■   C 

and 

^^C  =    -^F^^Es/^l^f^'^'-'^-^B/^Xl^J^^EB  +  ^Est'^^'-^^-^EB/^Xl^-^^^^F 

+  1   /v   [c,xp(-v  /V   )]AV 

CS   X2    '     CB   X2  ■■   CB 

Substituting  these  relations  into  (3.11)  for  the  two  branches 
yields 


AV   <'  -A I  H'   +AV   $  ~AI  H' 
EB  E    E  EB    CB  C    C  CB 


B 


AV   <J)  -{I   /V   [exp(-V   A   )]AV 
EB  E    ES   XI     '     EB   XI      EB 


-a  I   /V  [dxpi-V       A7      )]AV   }T 
R  CS   X2    '     CB   X2      CB   EB 


+AV   $  -{-a  I   /V   [exp(-V   /V   )]AV 

CB  C     F  ES   XI         EB   XI      EB 


+1   [exp(-V  /V   )-l]Aa 
ES         EB   XI        F 


+  1   /V   [exp(-V   /V   )1AV.  }"]! 
CS   X2         CB   X2      CB   CB 


Combining  like  terms  yields 


^V^Es/^Xlf^'^P^-^EB/^Xl^l^-B 
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F  ES   XI     '     EB   XI  -^  CB    E 


B 


+  {$^+a^I^^A^^[exp(-v^3A^^)]Y^3 


-I  /v  [(ixp(-v  A  )]^  }^v 

CS   X2         CB   X2     CB    CB 


-I^3[exp(-V^3A,,)-l]H',3Aa^._ 


(3.25) 


We  eliminate  the  terms  containing  AV    and  AV   by  defining 


EB 


CB 


the  adjoint  network  relationship  to  be 


=  Y 


¥ 


EB 


CB 

L.     J 


where 


Y  = 


i  A  e  xp  (-V  A  ) 

ES   XI    '     EB^  XI 


-"r^Cs/^X2^"^^^(-\b/^X2^ 


-^f^es/^xi^'^'^^-^eb/^xi^ 


I  A  exp(-v  A  ) 

CS'  X2   ^    CB^  X2 


which  are  the  expressions  describing  a  linear  two-port.   The 
y-parameter  matrix  of  this  two-port  is  the  transpose  of  the 
y-parameter  matrix  of  the  final  linearization  of  the  tran- 
sistor model  in  the  original  network.   We  may  transform  the 
remaining  term  of  (3.25)  to  be  dependent  on  A(l/Bp.)  instead 


of  Aa  by  realizing  that 


9a^/3(l/6p)  =  3(l+l/p^)  V9(l/Bp)  =  -a^, 


Therefore,  (3.25)  becomes 
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a^I^3[exp(-V^^A^i)-l]\B^(l/6p) 


(3.26) 


The  ac  transistor  models  are  linear  and  may  be  represented 
by  the  two-port  expression 


T 


7 


V 


1 1 


21 


12 


22 


1/ 


EB 


CB 


(3.27) 


We  will  now  perturb  the  i'-parameters  and  determine  an  ex- 
pression that  can  be  used  to  specify  the  adjoint  network 
representation  and  gradient  information  for  each  transistor 
model  individually.   Perturbing  the  V-parameters ,  subtracting 
(3.27),  and  neglecting  second-order  terms  yields 


AI  =    V      Al/   +1/   AV   +7   Al/   +1/   A^ 
E      11   EB   EB   11    12   CB   CB    12 


and 


AI      =    y      Al/      +1/      ^V  ^+V  ^h\l+\I  ^^V  ^^. 
C  21        EB       EB        21        22       CB        CB        22 


Substitution  of  these  expressions  into  (3.11)  yields 
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Al/   4>  -AI  H'   +AI/   'I'  -AI^T^^ 
EB  E    E  EB    CB  C    C  CB 


Al/   $  -fy'   Al/   +1/   AV 
EB  E    1 1   EB   EB   11 


+  y       Al/   +1/   AV   IH'   +AI/   4) 
12   CB   CB   la-"  EB    CB  C 


-W       Al/   +1/   AV  +y       Al/   +1/   AV^^^J^t'  ^ 
^21   EB   EB   21   22   CB   CB   22-"  CB 


which  can  be  rearranged  into  the  form 


^     E        11     EB        21     CB-"         EB 


*■     C        12     EB        22     CB-"        CB 


[1/       H*       A/       +1/       Y       A/ 
EB    EB        11        CB     EB        12 


+  1/       Y       AV       +1/       T       AV    ^]  . 
EB    CB        21        CB     CB        22-" 


(3.28) 


The  adjoint  network  representation  for  a  network  represented 
by  /-parameters  is  a  network  with  the  transpose  of  the 
original  V-parameter  matrix: 


•I- 


11        21 


12        22 


EB 


H' 


CB 


The  linearized  Ebers-Moll  model  can  be  described  by  the 
following  V'-parameters : 
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V        =  (3  +1)  [  (1/3  +jw/w   )a  I   /V   £xp(V  /V   ) 

11       F     '      F  -^  '^  Tl   F  ES   XI   '    Dl^  XI 


+jwC   /(1-V   A   )**(l/m  )], 
■^   01       Dl   01       ^  1  -■  ' 


i^   =  -a  I  /V  exp  (V  /V   )  , 

12       R  CS   X2   '    D2   X2 


V        =  -a  I  /V   exp (V  /V   ) , 

21       F  ES'^  XI   '    Dr   XI  ' 

and 


V       =  (3  +1)  [  (1/3  +jaj/w   )a  I   /V   exp(V  /V   ) 

22       R     "^   ^  R  -^  ^  T2   R  CS'  X2   '    02^  X2 


+jwC   /(1-V   /V   )**(l/m  )]. 
■^   02       D2'^  02     ^  /  2  -■ 

The  four  parameters  for  this  transistor  model  are  1/3  ,  l/co   , 
^  F      Tl 

V   ,  and  V^„.   An  expression  in  terms  of  perturbations  of 
Dl         D2  ^  ^ 

these  parameters  can  be  obtained  from  (3.28)  by  the  following 
conversion : 


^^1  "  [S^j^/3(l/3p)]A(l/3p)  +  [8V^^/8(l/aJ^^)]A(l/a)^^) 


+  [8S'   /SV   ]AV   +[9^^   /3V   ]AV   . 
11    Dl    Dl      11     D2    D2 

Expressions  for  A^   ,  AV   ,  and  A^    are  similar  and  yield 
^  12      21  22 
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-'EB^'EBt^'/V^^/\l^"p'Es/^Xl^^^^^D/^l^ 


-T   a  I   /V^  exp(V  /V   )  }AV  ^ 

CB  F  ES   XI       Dl   XI      Dl 


-V        {-H'   a  I   A^  £Xp(V   /V   ) 
CB    EB  R  CS   X2       D2   X2 


+M'   [  (1/a  +jw/a)   )a  I   A  .^^P  (V^^o/^vo) 
CB*-   '     R    -*  T2   R  CS   X2   '    D2   X2 


The  hybrid-pi  model  can  be  described  by  the  following 
/-parameters : 


V        =  -nq|i  |/(kT) , 

12  C 


V         -  -qll  I  (l+n)/(kT)  , 


and 


y        =  nqjl  |/(kT)(l/3  +ju)/w^^  +  l)H-ja3C 
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The  three  parameters  for  this  transistor  model  are  1/B  /  1/w   , 

F       Tl 

and  |l  |.   An  expression  in  terms  of  perturbations  of  these 
parameters  can  be  obtained  from  (3.28)  by  the 'following 
conversion: 

^^11  =  [9^i/8(l/Bp)]A(l/3p)  +  [3i^j/3(lA^j)]A(l/a3^^) 


-{[2)V     /d\i    MaIi 


Expressions  for  A/   ,  A^   ,  and  hV        are  similar  and  yield 
^  12      21  22 

-qll  |/(kT)(l/   T   -l-nl^   ^   )A(1/B  ) 
'  C  EB  EB    CB  CB        F 


-jwqll  |/(kT)(l/  "V      +T]V      4'   )A(l/a3   ) 

J  ^1  ^1/     ' '  EB  EB     CB  CB        Tl 


-q/(kT)[(l/3  +j(io/w  +l+n)l^  ^  -nl^  "i^ 

^i/  V   /  L  V  /  Hp  J  /  ^^     "  EB  EB    CB  I 


EB 


-(n+Df     T     +n(l/3  +jaj/w     +1)1/     4'     ]a|i   ].  (3.30) 

EB  CB     '  F  -^    Tl      CB  CB      c' 

The  only  other  elements  that  are  present  in  the  original 
network  are  independent  voltage  and  current  sources,  some  of 
which  are  zero-valued  and  introduced  solely  for  the  purpose 
of  measuring  a  response.   An  independent  current  source  is 
required  for  every  desired  voltage  response  and  an  indepen- 
dent voltage  source  is  required  for  every  desired  current 
response  because  they  will  be  needed  as  sources  in  the 
adjoint  network  even  if  they  are  zero-valued  in  the  original 
network.   We  assume  that  the  values  of  the  independent  sources 
remain  constant  so  that 
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Expression  (3.11)  for  ac  considerations  reduces  to 

-IAI,,vI^^^  =  - 1  ^jA^j  (3.31) 

V  I  J=l 

where  n  is  the  number  of  designable  parameters,  r   is  the  jth 
designable  parameter,  and  G      is  the  corresponding  coefficient 
which  can  be  found  in  (3.12)  through  (3.18),  (3.20),  (3.29), 
and  (3.30).   For  dc  considerations,  (3.11)  reduces  to 

V  I  J=l 

where  G   can  be  found  in  (3.12),  (3.13),  (3.18),  (3.19),  and 
(3.26).   We  can  associate  independent  sources  in  the  original 
network  v/ith  independent  sources  in  the  adjoint  network  of 
the  same  type.   The  values  of  these  sources  are  still  to  be 
found. 

We  nov7  consider  the  ac  performance  function  (3.6). 
Consideration  of  other  ac  performance  functions  yields  sim- 
ilar results.   The  change  in  error  due  to  the  change  in 
designable  elements  is,  to  first  order, 

AE,_  =  I    j        RedW  (u)  ([1/  (u)  |/|P  (u)  |-1)]P  ^ 

AC     J    ^  I  I  I 

I 

•W^  (u)l'^  (u)/[|l/^  (u)  I  |l/^  (u)  |]Al/  (u)}du 


+  1  /    Re{[W^^(u)  (|l^,(u)  |/|I^,(u)  |-1)]P-1 

I 
I 


V   u         V      '  V     "  '  V 


* 


•W^(u)I^(u)/[|l  (u)||l  (u)|]AT  (u)}du      (3.33) 
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where  *  indicates  complex  conjugate.   Comparison  of  (3.31) 
and  (3.33)  reveals  that  if  all  current  sources  in  the  ad- 
joint network  corresponding  to  desired  responses  are  set  to 

$j(u)  =  [Wj(u)  (|l/j(u)  |/ll/j(u)  |-1)]P~1 

.[W^  (u)l/*(u)/(|U^(u)  I  Il/^(u)  I)], 

all  voltage  sources  in  the  adjoint  network  corresponding  to 
desired  responses  are  set  to 

'i'^(u)  =  -[W^(u)  (Il^(u)  |/|l^(u)  |-1)]P-1 

•[W^(u)T*(u)/(lT^(u)|l7^(u)|)], 

and  all  remaining  sources  in  the  adjoint  network  are  set  to 
zero,  then 

^^AC  =   I   ^/    Re(-Gj)du}Arj. 
J=l    Uj 


Since 


J=l 


the  desired  gradient  components  result: 

"f 

9E   /8r   =  /    Re(-G  )du. 
AC    J    •'  „        J 


u 
I 


A  similar  process  specifies  the  adjoint  network  excitations 
for  do  considerations  to  be 
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4>j(T)  =  [Wj(T)  (Vj  (T)/Vj(T)-l)]P  ^Wj(T)/Vj(T) 


and   ^ 


^     (T)  =  -[W^  (T)  (I^  (T)/£  (T)-1)]P  ^  W  (T)/i  (T) 
V  V       V       V  V       v 

when  using  performance  function  (3.3).   The  dc  gradient 
components  are 

I 
Since  the  total  error  is  the  sum  of  the  ac  error  and  - 

the  dc  error,  the  total  gradients  are  the  sum  of  the  ac 

gradients  and  dc  gradients: 

We  have  determined  first-order  partial  derivatives  of 
the  performance  function  with  respect  to  each  netv;ork  param- 
eter.  Since  each  of  these  parameters  has  a  temperature 
dependence  of  the  form 

r   =  r ' [1+a  (T-T  ) ] 

where  r   is  the  component  of  the  parameter  that  must  be 
varied,  we  determine  gradient  information  v;ith  respect  to  r 
and  transform  it: 


8E/3r   =  [8E/9r  ] [3r  /dr    ] 


=  [l+aj(T-T^)][3E/9rj]. 
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It  has  been  shown  that  the  gradients  of  all  network 
paraineters  of  interest  are  appropriately  weighted  products 
of  voltages  and  currents  of  the  original  and  adjoint  networks. 
All  gradient  information  can  be  determined  upon  performing 
two  ac  analyses  per  frequency  and  two  dc  analyses  per  temper- 
ature.  (It  should  be  noted  that  the  adjoint  network  is 
linear  and  requires  only  one  linear  dc  analysis.)   The  gradi- 
ents obtained  by  this  method  have  been  shown  to  be  exact.® 

Accurate  Integration  of  the  Performance  Function  and  Gradients 

The  evaluation  of  the  performance  functions  and  gradients 
discussed  above  requires  integration  of  functions  of  network 
responses.   Because  of  the  numerical  nature  of  the  network 
functions,  the  ac  performance  function  and  the  ac  gradient 
components  must  be  numerically  integrated  over  the  frequency 
range  and  the  dc  performance  function  and  dc  gradient  compo- 
nents must  be  numerically  integrated  over  the  temperature 
range.   A  set  of  frequencies  and  a  set  of  temperatures  must 
be  selected  by  the  designer  or  the  design  algorithm  that 
accurately  represent  the  entire  domain  of  operation.   The 
quality  of  the  final  network  depends  on  the  accuracy  of  the 
numerical  integration  algorithms. 

All  of  the  functions  to  be  integrated  are  explicit 
functions  of  frequency  or  temperature  and  it  is  not  possible 
for  the  numerical  integral  to  be  unstable  unless  the  function 
itself  is  unstable. 
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The  dc  performance  function  and  do  gradient  components 
are  generally  smooth  functions  of  temperature  that  can  be 
integrated  accurately  by  selecting  three  or  more  fixed  tem- 
peratures and  using  trapezoidal  rule.   Desired  dc  responses 
are  usually  constant  voltage  and  current  specifications. 
Actual  dc  responses  are  usually  linear  within  the  range  of 
temperatures  over  which  the  element  models  are  valid. 

The  ac  performance  function  and  ac  gradient  components 
may  be  rapidly  varying  functions  of  frequency  and  in  such 
cases,  a  fixed  set  of  frequencies  will  not  be  adequate  unless 
a  large  number  of  frequencies  are  specified..  An  algorithm 
that  automatically  selects  frequencies  in  order  to  maintain 
a  specified  accuracy  enables  the  integration  to  be  performed 
with  a  smaller  num.ber  of  frequencies.   The  network  response 
shown  in  Fig.  3-3  is  typical  of  the  type  of  response  that 
best  utilizes  the  flexibility  of  an  automatic  frequency 
selection  algorithm.   The  response  from  03   to  co   is  unchanged 

A       B 

and  few  frequencies  are  necessary  to  numerically  integrate 
this  portion  of  the  response  accurately.   On  the  other  hand, 
the  response  from  w   to  to   varies  considerably  and  many 

B       C 

frequencies  must  be  used  to  accurately  integrate  this  portion 
of  the  response. 

Integration  of  the  dc  functions  is  accomplished  by  use 
of  trapezoidal  rule  at  a  set  of  tem.peratures  specified  by  the 
user.   The  formula  used  is 
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OUT 


Fig-  3-3 


Network  response  over  a  frequency  range 
with  peaking  where  solid  line  is  actual  response 
and  dashed  line  is  desired  response 
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T 

N 

/    f(T)dT  =  f  (T^)  (T2-T^)/2  +  f  (T^)  (T^-T  ^  ) /2+ .  .  . 
1 

+  f(T    )  (T  -T    )/2+f(T  )  (T  -T    ) /2     (3.34) 
N-1     N   N-2^       N    N   N-1 

where  N  is  the  number  of  temperatures  and  the  function  f(T) 
can  be  either  the  performance  function  or  a  gradient  compo- 
nent.  Accuracy  may  be  increased  by  specifying  more  tempera- 
tures, but  there  is  seldom  need  for  more  than  ten. 

Three  different  numerical  integration  algorithms  were 
tested  for  evaluation  of  the  ac  performance  function  and 
gradients.   Two  of  the  algorithms  were  found  to  be  more 
efficient  for  the  functions  of  interest  here. 

The  most  straightforward  algorithm  uses  trapezoidal 
rule  at  a  set  of  frequencies  specified  by  the  user.   As  with 
the  dc  case,  (3.34)  is  used  with  angular  frequency  as  the 
variable  of  integration  instead  of  temperature.   It  was 
found  that,  except  for  some  simple  networks,  a  large  number 
of  frequencies  had  to  be  specified  to  achieve  acceptable 
results. 

The  second  integration  algorithm  also  employed  trape- 
zoidal rule,  but  included  an  automatic  frequency  selection 
feature.   Tlie  frequency  selection  process  is  based  on  a 
halving  or  doubling  of  the  step  size  depending  on  the  degree 
to  which  a  specified  error  criterion  is  satisfied.   The 
error  criterion  compares  the  partial  integral  evaluated  over 
a  double  step  size  to  the  sum  of  two  partial  integrals,  each 
evaluated  over  a  single  step  size,  and  indicates  whether  the 
step  size  is  suitable. 
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To  begin  the  integration,  three  analyses  are  performed. 
The  first  analysis  is  performed  at  the  initial  frequency  of 
the  range  and  all  function  values  for  that  frequency  are 
saved.   The  frequency  is  updated  by  adding  to  it  one- tenth 
of  the  entire  frequency  range.   (When  the  frequency  is 
scaled  logarithmically,  the  step  size  is  logarithmic  also.) 
One-tenth  of  the  range  is  specified  to  be  the  initial  step 
size  and  the  maximiim  step  size  because  it  is  small  enough  to 
prevent  most  aberrations  from  being  overlooked  and  yet  it  is 
large  enough  to  prevent  the  use  of  too  many  frequencies  when 
they  are  unnecessary.   After  analysis  at  the.  second  frequency, 
all  function  values  are  stored  and  a  third  frequency  is 
obtained  by  adding  the  step  size  to  the  previous  frequency. 
Analysis  at  this  frequency  allows  two  values  of  the  integral 
to  be  calculated  and  compared: 

Q  =  |h/2[f(a)  )+2f(a)  )+f(a3  )]-h[f(a;  )+f(w  )  ]  | /A   (3.35) 

1  ^  ^  J.  <3 

where  h  is  the  step  size, 

(jo   is  the  first  frequency 

w   is  the  second  frequency, 

0)   is  the  third  frequency, 

A  is  the  smallest  maximum  value  previously  assigned  to 
the  integral  of  f(w),  and 

Q  is  the  measure  of  the  acceptability  of  the  current 
step  size. 
If  A  is  zero,  as  it  would  be  on  the  first  comparison  of  the 
first  analysis,  the  normalization  of  the  difference  is  not 
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performed.   Since  many  functions  are  being  integrated  simul- 
taneously, the  values  of  all  Q's  are  averaged  to  obtain  Q^ 
which  is  used  to  judge  the  acceptability  of  the  step  size. 
If  Q  >2e  where  e  is  a  specified  error  criterion,  the  step 
size  is  unacceptably  large  and  is  halved  unless  it  is  already 
at  the  minimum  step  size.   If  £;lQ^l2e,  the  step  size  is 
acceptable  and  its  value  is  retained.   If  Qf^<^'    the  step  size 
is  unacceptably  small  and  is  doubled  unless  it  is  already  at 
the  maximum  step  size.   When  the  step  size  is  retained  or 
increased,  integration  proceeds  from  the  third  frequency;  if 
the  step  size  is  decreased,  integration  resumes  at  the  first 

frequency.     •   ~ 

The  ac  performance  function  and  gradients  are  explicit 
functions  of  frequency.   By  explicit  functions,  we  mean  that 
the  value  of  such  a  function  at  a  specific  frequency  is  not 
dependent  on  the  value  of  the  function  at  another  frequency. 
Thus,  if  a  slight  inaccuracy  is  allowed  early  in  an  integra- 
tion, this  inaccuracy  will  not  affect  the  remainder  of  the 
integration.   Since  only  explicit  functions,  are  considered, 
error  accumulation  is  not  examined  in  the  halving  and  dou- 
bling integration  scheme  and  a  straightforward  algorithm 
results . 

The  third  integration  algorithm  is  a  direct  descendant 
of  the  variable  step  size,  variable  order  algorithm  developed 
by  Gear.^^   This  algorithm  is  a  predictor-corrector  multi- 
step  numerical  integration  scheme  based  on  tv/o  variable-order 
equations.   The  first  equation  is  used  to  predict  the  partial 
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integral  for  a  specified  step  size  and  order  and  the  second 
equation  is  used  iteratively  to  correct  the  step  size  and 
order  to  reduce  the  error  in  integration  be.lovj  a  specified 
value.   Gear's  algorithm  is  intended  to  integrate  implicit 
functions  and  has  an  elaborate  error  control  mechanism  built 
into  it.   Because  of  the  complexity  of  this  algorithm,  it  is 
the  most  inefficient  of  the  three  algorithms  although  it 
generally  performs  acceptably.   An  error  criterion,  e    ,  is 
specified  for  this  algorithm  which  is  similar  to  the  error 
criterion  of  the  other  variable  step  size  algorithm  but  not 
identical. 

A  problem  common  to  both  variable  step  size  algorithms 
is  the  possibility  of  obtaining  a  slightly  different  value  of 
the  performance  function  on  two  passes  even  though  the  net- 
work parameters  are  the  same.   This  inconsistency  occurs 
because  the  value  of  A  in  (3.35)  is  constantly  being  updated. 
Since  the  optimization  algorithm  becomes  confused  at  this 
fluctuation,  a  reasonably  tight  error  criteria  should  be 
maintained  to  reduce  the  significance  of  this  possible  varia- 
tion. 

The  three  algorithms  were  compared  through  some  typical 
design  statistics  obtained  through  use  of  the  computer  pro- 
gram FROLIC.   The  first  network  to  be  considered  is  a  seven- 
pole,  LC,  passive,  low-pass  filter  with  initial  parameter 
values  as  shown  in  Fig.  3-4 (a).   The  desired  response  and 
the  initial  actual  response  to  a  constant  one  ampere  excit- 
ation are  shown  in  Fig.  3-4 (b) .   The  variable  network  para- 
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(a)  Seven-pole  passive  lov7-pass  filter 

(b)  Desired  and  initial  actual  responses 
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meters  are  the  four  capacitances  and  the  three  inductances. 
For  all  of  the  designs,  performance  function  (3.6)  was  used 
with  p=2 ,  optimization  was  by  means  of  the  conjugate  gradients 
algorithm,  the  parameters  were  logarithmically  scaled,  and 
the  parameter  constraints  were  set  so  that  they  did  not  affect 
the  design  process.   The  values  of  e  and  e   are  varied  for 
the  variable  step  size  algorithms  and  the  nun^oer  of  frequen- 
cies is  varied  for  the  fixed  frequencies  algorithm.   Compar- 
isons are  made  in  Table  3-1  of  average  time  per  analysis, 
average  time  per  frequency,  average  number  of  frequencies 
per  analysis,  and  whether  the  desired  result  was  obtained. 
The  fixed  frequencies  were  chosen  so  that  they  were  more 
dense  in  the  vicinity  of  the  break  frequency.   It  is  shown 
that  more  frequencies  m.ay  be  used  in  the  fixed  frequency 
algorithm  because  it  has  no  overhead  for  frequency  selection; 
although  fewer  frequencies  may  be  just  as  acceptable  because 
of  their  appropriate  placement  with  the  variable  step  size 
algorithm.   It  is  also  shown  that  the  trapezoidal  rule 
variable  step  size  algorithm  is  more  efficient  than  the  Gear 
algorithm  for  this  type  of  problem. 

For  a  second  example,  we  design  the  response  of  a  series- 
series  triple  feedback  amplifier  with  parameters  and  initial 
responses  shown  in  Fig.  3-5.   Three  designs  are  run  with  each 
of  the  three  integration  algorithms.   For  all  of  tlie  designs, 
performance  function  (3.6)  was  used  with  p=2,  optimization 
was  supervised  by  the  conjugate  gradients  algorithm  and  the 
parameters  were  logarithmically  scaled.   The  values  of  e  and 
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Fig.  3-5 
Series-series  triple  feedback  amplifier 
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c   are  set  to  ,01  except  in  one  run  v/here  Gear's  algorithm 
was  being  used  and  unacceptable  results  were  obtained.   In 
this  case,  a  value  of  e  =.005  was  used  and  satisfactory 
results  were  obtained.   The  value  .01  has  been  found  to  be 
acceptable  for  most  problems.   Comparisons  are  made  in 
Table  3-2  of  average  time  per  analysis,  average  time  per 
frequency,  average  number  of  frequencies  per  analysis,  and 
whether  the  desired  result  was  obtained.   The  first  desired 
response  was  a  flat  gain  of  fifty  over  a  bandwidth  of  one 
hertz  to  a  hundred  megahertz.   All  capacitors  and  resistors 
in  the  network  except  the  source  resistor  were  allowed  to 
vary.   The  results  shov;ed  all  algorithms  to  be  capable  of 
simultanviously  integrating  the  seventeen  functions,  but  the 
variable  step  size,  trapezoidal  rule  algorithm  is  superior. 
If  eight  of  the  capacitors  are  removed  from  the  netv;ork  and 
only  C  ,  C  ,  and  the  resistors  remain  as  network  parameters, 
similar  results  are  obtained.   If  all  capacitors  are  removed 
from  the  network  and  the  desired  bandv;idth  reduced  slightly, 
the  two  trapezoidal  rule  algorithms  are  successful.   Gear's 
algorithm  had  to  have  a  reduced  value  of  e '  to  achieve  the 
desired  result. 

The  fixed  frequency  approach  is  a  desirable  algorithm 
for  some  problems  and  v.'ill  work  for  all  problems  if  enough 
frequencies  are  specified.   It  has  the  advantage  of  always 
yielding  the  same  result  for  identical  functions.   The 
variable  step  size  trapezoidal  rule  algorithm  has  performed 
well  and  frees  the  user  of  the  chore  of  predicting  areas  of 
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the  frequency  range  where  problems  might  occur.   This  algo- 
rithm is  superior  to  Gear's  algorithm  because  it  uses  fewer 
frequencies  and  is  more  efficient  when  integrating  the  explic- 
it functions  considered  here. 

An  additional  feature  of  a  variable  step  size  algorithm 
is  that  if  aberrations  in  the  frequency  response  move  around 
within  the  frequency  range,  the  algorithm  can  adjust  to 
their  movement.   If  these  aberrations  disappear  during  the 
course  of  the  design,  the  algorithm  can  adjust  for  this  also. 


CHAPTER  4 
TOPOLOGY  DESIGN  THROUGH  ELEMENT  GROWTH  AND  REMOVAL 

Parameter  adjustment  to  obtain  desired  network  responses 
solves  only  part  of  the  network  design  problem.   An  initial 
network  configuration  must  be  determined  that  has  the  possi- 
bility of  yielding  the  desired  responses-   In  this  chapter, 
a  practical  approach  to  the  design  of  network  topology  is 
discussed  with  emphasis  on  the  growth  and  removal  of  tran- 
sistors to  satisfy  ac  specifications. 

Resistor  and  Capacitor  Growth  and  Removal 

A  fixed  network  configuration  limits  the  flexibility  of 
an  automated  design  procedure.   The  ability  to  add  and  remove 
elements  allows  nev7  configurations  to  be  developed  and  old 
configurations  to  be  improved.   The  growth  and  removal  of 
linear  elements  has  been  discussed  by  Director  and  Rohrer 
in  conjunction  with  the  adjoint  network  method  of  gradient 
calculation. 

The  gradient  for  a  resistor, 

8E/8R  =  /    Re(-I  $  )  [1+a  (T-T  )  ] dto 
'  '  R  R       R      0 

+/    -^r'^r^^+^r^'^-'^O^^^'^' 

does  not  explicitly  depend  on  the  value  of  resistance.   If 
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■R-n   T    ^  .  I  ,  and  (b      can  be  calculated  and  the  desirability 
of  "growing"  a  resistor  in  a  short  circuit  is  determined. 
The  gradient  for  a  conductor,  actually  a  resistor  with  G=l/R, 
is 

8E/3G  =  /    Re(l/g7^)  [l+a^(T-TQ)  Jdu 

+  /    V^ij;^[l+a^(T-T^)]dT 

which  is  not  dependent  on  the  value  of  conductance.   If  G=0, 
1/    vi'(-,  V^,    and  '^^    can  be  calculated  and  the  desirability  of 
growing  a  conductor  in  an  open  circuit  is  determined.   Although 
a  resi-tnr  and  a  conductor  are  simple  inverses  of  each  other 
at  any  non-zero  values  of  resistance  and  conductance,  it  is 
the  zero  values  that  are  of  interest  when  considering  growth. 
Thus  two  representations  of  a  resistive  branch  are  considered, 
one  to  grow  from  a  short,  circuit  and  one  to  grow  from  an 
open  circuit.   Similarly,  the  gradient  for  a  capacitor  will 
indicate  possible  growth  of  a  capacitor  in  an  open  circuit 
and  the  gradient  for  an  elastance  indicates  possible  growth 
of  a  capacitor  in  a  short  circuit. 

With  indicators  available  to  determine  if  a  resistor  or 
a  capacitor  should  be  added  in  any  lead  or  between  any  pair 
of  nodes,  it  initially  appears  to  be  possible  to  check  all 
locations  for  resistor  and  capacitor  growth.   Actual  appli- 
cation of  this  procedure  would  require  more  time  than  is 
feasible  and  an  alternate  approach  is  used:   possible  growth 
locations  are  specified  by  entering  components  with  minimum 
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parameter  values.   It  is  not  possible  to  have  components 
with  zero-valued  parameters  because  of  the  general  branch 
model  assumed  earlier. 

Since  zero-valued  parameters  are  not  used,  only  resist- 
ances and  capacitances  need  to  be  considered.   A  minimum- 
valued  conductance  is  equivalent  to  a  maximum- valued  resist- 
ance, and  similarly,  a  minimurri-valued  elastance  is  equivalent 
to  a  maximum- valued  capacitance. 

Resistor  and  capacitor  removal  is  presently  simulated  by 
the  adjustment  of  their  parameters  to  their  lower  or  upper 
limits.   VJhen  these  parameters  encounter  a  boundary,  the 
elements  are  not  removed  from  consideration  as  they  may  again 
become  necessary  to  achieve  the  desired  responses  and  "grow" 
by  their  parameters  moving  av/ay  from  the  boundaries.   This 
procedure  has  proven  to  be  satisfactory. 

Diode  Growth  and  Removal 

The  gradient  for  a  diode  (the  contribution  due  to  the 
junction  voltage  is  omitted  from  growth  consideration)  is 

03 

F 

+/    [(ixp(V^A^)-l]il'j5[l+ap(T-T^)]dT 

where  1/  ,  H*  ,  V  ,  and  \b      are  voltages  across  the  junction  and 
D    D    D        D 

the  expression  in  the  second  integrand  is  a  simplified  form 
of  (3.19),   Since  the  expression  contains  no  dependence  on 
the  parameter  I  ,  it  is  possible  to  grow  a  diode  between  any 
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tv70  nodes.   The  problem  with  growing  a  diode  is  that  the 
gradient  is  actually  indicating  conductor  growth.   Thus  a 
diode  is  grown  as  a  variable  conductance  and  its  voltage 
setting  characteristics  which  are  a  consequence  of  the  diode's 
nonlinearity  are  ignored.   An  alternate  indicator  to  measure 
the  desirability  of  a  fixed  voltage  has  not  been  determined 
so  diode  growth  and  removal  are  not  considered. 

Transistor  Growth  and  Removal 


Bipolar  junction  transistors  have  become  the  most  com- 
monly used  components  in  integrated  circuits  because  tran- 
sistors are  small,  dissipate  little  power,  and  pex'form  many 
functions.   Transistor  growth  and  removal  is  therefore  a 
relevant  subject  to  consider.   In  this  work,  dc  considerations 
are  not  used  to  determine  transistor  growth  or  removal.   The 
hybrid-pi  transistor  model  is  used  because  it  is  convenient 
and  yields  results  comparable  to  those  obtained  through  use 
of  other  models. 

Initially,  we  approach  transistor  growth  in  a  manner 
similar  to  that  applied  to  other  components:   we  investigate 
the  gradient  expressions  for  a  zero-valued  transistor.   A 
zero-valued  transistor  may  be  envisioned  as  one  with 

1/Bp  =  1/co^j  =  \l^\    =    0.  ■ 

In  this  case,  the  only  ac  transistor  gradient  component  not 
identically  equal  to  zero  is 

8E   /8|l  i  =  q/(kT)/    Re[V      <}/   -1/   y   ]dw 
AC     C  \^  EB  EB   EB  CB 

I 
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V7hich  can  be  rewritten   as 

9E^^/8|I^|    =q/(kT)/         Rell/^^^F^^Jdu)  (4.1) 

where  n=0  for  a  nonexistant  transistor.   The  gradient  of  a 
voltage-controlled  current  source  shown  in  Fig.  4-1  is  known 
to  be 

0) 

F 

3E   /Sg   =  /    Re[l/   H'   ]  du . 

AC    M      ,,        EB  EC 

a3j 

Since  g  =[q/ (kT) ]  I  I  I  ,  the  gradient  computation  associated 
with  a  zero-valued  transistor  is  the  same  as  that  for  a 
voltage-controlled  current  source. 

At  this  juncture  it  is  important  to  assess  the  practi- 
cality of  v;hat  has  been  indicated  in  the  previous  paragraph. 
The  desirability  of  growing  a  voltage-controlled  current 
source  between  three  distinct  nodes  is  easily  determined. 

Since  the  voltages  V        and  Y    can  always  be  ascertained, 

^     EB       EC        . 

even  for  the  case  of  9f.=0 ,    the  gradient  expression  (4.1)  may 
be  evaluated.   But  for  many  situations  it  proves  more  con- 
venient to  consider  the  possibility  of  growing  transistors 
simultaneously  in  a  short  circuit  and  across  an  open  circuit. 
For  instance,  consider  the  single  common-emitter  stage  of 
Fig.  4-2 (a).   An  additional  common- emitter  stage  might  be 
required  if  more  gain  is  needed.   If  so,  another  common- 
emitter  stage  should  then  be  grown  as  indicated  in  Fig.  4-2  (b) . 

In  general,  there  are  three  possibilities  for  transistor 
growth  as  illustrated  in  Fig.  4-3.   Observe  that  for  the 
common- emitter  configuration,  the  short  circuit  is  between  the 
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Fig.  4-1 

Voltage-controlled  current  source  configuration 
that  is  equivalent  to  a  zero-valued  transistor 
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Fig.  4-2 

(a)  Single-stage  transistor  amplifier  with  growth  site 

(b)  The  growth  of  an  additional  conimon-emitter  stage 
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Fig.  4-3 
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Transistors  grown  siraultaneously  in  a  short  circuit 
(between  terminals  1  and  2)  and  across  an  open  circuit 
(between  terminals  1  or  2  and  3) .   The  three  possible 
ways  in  which  a  transistor  may  be  grown  are  shown. 
Common-emitter  stage. 

(c)  Common-base  stage. 

(d)  Comnion-col lector    staae. 
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base  and  collector  nodes  and  the  gradient  component  (4.1) 
reduces  to  that  of  a  conductance  since  I'  =\l        and  ^      =H' 

EC   EB        EB   EC 

Moreover,  the  gradients  associated  v/ith  the  common-base  and 

common-collector  are  identically  zero  since  I'   =¥   =0  and 

^  EC   EC 

f  ^g='l''^g=0 ,  respectively.   Thus  the  choice  of  a  voltage- 
controlled  current  source  to  indicate  transistor  growth  is  a 
poor  one.   We  are  therefore  forced  to  alter  our  approach. 

It  is  reasonable  to  ask  what  causes  a  need  for  transistors, 
the  answer  being  the  desirability  for  voltage  gain  and/or 
current  gain.   It  can  be  shown  that  in  an  ideal  situation, 
each  transistor  configuration  yields  gain.   The  gain  for  each 
stage  can  be  determined  by  use  of  an  ideal  mid-frequency 
transistor  model  in  each  of  the  three  possible  configurations 
as  shown  in  Fig.  4-4.   The  voltage  and  current  gain  expressions 
for  each  configuration  are  shovm  in  Table  4-1.   The  common- 
emitter  configuration  yields  a  large  voltage  gain  for  typical 
values  of  transconductance  and  load  resistance  and  a  large 
current  gain  for  typical  values  of  beta.   Both  the  voltage 
and  current  gains  of  the  common-emitter  include  phase  inver- 
sion that  can  be  ignored  for  a  magnitude  fit  (3.5)  but  must 
be  considered  for  a  complete  fit  (3.2).   The  common-base 
configuration  also  yields  a  large  voltage  gain  but  it  yields 
no  current  gain.   The  common-collector  configuration  yields 
no  voltage  gain  and  large  current  gain.   Therefore,  transistor 
growth  should  be  dependent  on  the  need  for  voltage  and/or 
current  gain. 
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Fig.  4-4 

Networks  used  to  determine  voltage  and  current  gain 
for  each  transistor  configuration 

(a)  Coniraon-emitter 

(b)  Common-base 

(c)  Common-collector 
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Gradients  that  indicate  growth  of  controlled  sources 
also  indicate  need  for  voltage  and  current  gain.   A  voltage- 
controlled  voltage  source  connected  as  shown  in  Fig,  4-5 (a) 
simulates  the  voltage  gain  of  the  common-base  configuration 
while  a  current-controlled  current  source  connected  as 
shown  in  Fig,  4-5 (b)  simulates  the  current  gain  of  the 
common-collector  configuration.   Both  a  voltage-controlled 
voltage  source  and  a  current-controlled  current  source 
connected  as  shown  in  Fig.  4-5(c)-{d)  simulate  the  voltage 
and  current  gains  of  the  common-emitter  configuration  for 
magnitude  fit  as  in  (3,5)  and  complete  fit  as  in  (3.2), 
respectively.   The  gradient  components  associated  v/ith  each 

of  these  configurations  are'' 

Dp 
8E   /8y  =  /    Re[l/  $  Ido) 

for  the  voltage  gain, 


"i 


to 

F 

9E^^/96  =  /    Re[-I^H'Jda3 


for  the  current  gain, 


dE      /8(y3)  E  dE      /8y+8E   /Bg 
AC  AC        AC 


Wp 


=  /    {Re[l/^<I'^]+Re[-I^H' J  }da) 
for  voltage  and  current  gain  for  a  magnitude  fit,  and 


Wp 


SE   /3(yB)  =  /    {Re[-l/,$  J+Re[T  ^    ]  }da) 

I 


Tl 

+ 

I/. 


^cv  W^cv=0 


yl^vcv 
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l2  = 
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i/2=(y+i)i^i 
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1^101=0 
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1/2=1/, 


(b) 


7p=-(3+l)7 


l'2=(y  +  l)l/j 


(c) 


^1   Ti(-j 


+   1/^^=0 


^cv=0 


<}- 


l2=-(l-3)Ti 


J)   VvW^^ICI     ^2=(-^-^)^l 


(d) 

Fig.  4-5 

Controlled  source  configurations 

(a)  Voltage  gain 

(b)  Current  gain 

(c)  Both  voltage  and  current  gain  with  inversion 
for  magnitude  fit 

(d)  Both  voltage  and  current  gain  with  inversion 
for  coiTiplete  fit 
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for  voltage  and  current  gain  for  a  complete  fit.   The  gradients 
for  corPLmon-base  and  common-collector  are  not  changed  by  use 
of  magnitude  or  com^plete  fit  because  the  desired  response  is 
positive  in  both  cases.   For  a  common-emitter,  the  desired 
response  must  change  signs  to  change  from  one  to  the  other  so 
the  gradient  expressions  must  also  differ  in  sign.   Observe 
that  for  the  zero-valued  case  (y=B=0) ,  each  of  these  gradients 
exist  and,  moreover,  each  of  the  models  of  Fig.  4-5  reduce  to 
the  appropriate  initial  model  of  Fig.  4-3(a).   Furthermore, 
after  two  circuit  analyses,  one  on  the  original  netv7ork  and 
the  other  on  the  adjoint  network,  it  is  not  only  possible  to 
determine  if  a  transistor  should  be  grown,  but  the  appropriate 
configuration  is  also  indicated  as  shown  in  Tables  4-2  and 
4-3. 

If  transistor  growth  is  indicated,  appropriate  initial 
param.eter  values  must  be  defined  to  simulate  as  closely  as 
possible  the  short  and  open  relationship  previously  in  effect. 
The  dc  bias  parameters  are  directly  related  to  the  elements 
connecting  each  pair  of  nodes.   Therefore,  appropriate  values 
of  these  parameters  may  not  be  at  boundaries  so  that  a  con- 
dition of  all  shorts  or  all  opens  is  avoided;  the  initial 
values  for  these  parameters  are  chosen  to  be  the  values 
specified  as  input  by  the  user  for  all  three  configurations. 

The  initial  values  of  l/[3   and  l/a>    for  the  comm.on-col lector 

F         Tl 

configuration  are  chosen  to  be  the  upper  limits  for  these 
parameters  in  order  to  simulate  the  short  from  base-to-emitter. 
The  value  of  n  is  typically  very  small  so  the  base-to-collector 
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Table  4-2 


Transistor  Growth  Indicators  for  Magnitude  Fit 


8E^^/3y 

8E^^/33 

Transistor  Configuration 

>0 

>0 

No  Transistor  Needed 

<0 

>0 

Grow  Common-Base  Stage 

>  0 

<0 

Grow  Common-Col] ectcr  Stage 

<0 

<0 

Grow  Common-Emitter  Stage 

JLIO 


Table  4-3 


Transistor  Growth  Indicators  for  Complete  Fit 


3E^^/3y 

3E^^/8e 

Transistor   Configuration 

>0 

>0 

Grow  Cornmon-Emitter   Stage 

<0 

>0 

Grow   CoiTunon-Base   Stage 

^  >o 

<0 

Grov;   Common-Collector   Stage 

<0 

<0 

No  Transistor  Needed 

Ill 


connection  does  not  simultaneously  become  a  short  circuit. 
For  the  other  two  configurations,  initial  values  of  1/3^  and 
1/co    are  chosen  to  be  the  lower  limits  in  order  to  simulate 
the  necessary  open  circuits;  the  short  circuit  is  difficult 
to  simulate  in  these  cases. 

The  ability  to  add  transistors  to  a  network  must  be 
complemented  by  the  ability  to  remove  transistors  that  are 
no  longer  desirable.   It  is  not  sufficient  to  allow  transistor 
parameters  to  remain  at  a  boundary  as  such  a  transistor  still 
has  significant  influence  on  the  network  response.   A  simple 
approach  to  transistor  removal  is  now  suggested. 

The  transistor  is  no  longer  useful  for  satisfaction  of 
ac  specifications  when  the  voltage  and/or  current  gain  it 
provides  is  not  desired.   An  indication  of  this  condition  can 
be  obtained  by  the  dererraination  of  the  desirability  of  con- 
ductor growth  between  the  terminals  of  a  transistor.   Gradients 
for  the  zero-valued  conductors  shown  in  Fig.  4-6  can  be  cal- 
culated and  used  as  removal  indicators.         ~ 

When  a  transistor  has  been  removed,  it  must  be  replaced 
by  the  short-open  combination  from  which  transistors  grow. 
The  location  of  the  short  must  be  determined  from  the  removal 
indicators.   This  can  be  done  as  follov;s.   First  removal 
requires  that  all  three  of  the  gradients  for  the  conductors 
in  Fig.  4-6  be  negative  (indicate  conductor  growth),  and 
that  the  conductor  with  the  most  negative  gradient  be  re- 
placed by  the  short  circuit.   This  method  is  simple  and  yet 
general.   If  the  rem.oval  of  a  transistor  causes  the  value  of 
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the  performance  function  to  increase,  the  transistor  is  re- 
stored and  other  transistors  are  checked  for  removal. 

Transistor  growth  indicators  may  be  determined  for  all 
possible  combinations  of  a  short  and  an  open.   This  approach 
would  consume  large  amounts  of  time  and  is  therefore  dis- 
carded in  favor  of  requiring  possible  growth  sites  to  be 
specified  a  priori.   Upon  transistor  removal,  the  location 
of  the  removed  transistor  becomes  a  growth  site.   Similarly, 
all  transistors  that  are  grown  automatically  become  removable. 

After  each  growth  or  removal,  a  design  is  performed  and 
a  minimum  is  obtained.   If  the  value  of  the  ac  portion  of  the 
performance  function  at  the  end  of  the  design  is  larger  than 
it  was  before  the  action  was  taken  on  the  transistor,  the 
action  is  reversed  and  the  old  parameter  values  are  restored. 

The  desirability  of  transistor  growth  is  checked  before 
transistor  removal  is  considered  because  it  is  more  often 
necessary  to  supply  missing  gain  than  it  is  to  remove  un- 
desirable gain.   It  is  not  possible  to  compare  growth  indi- 
cators with  removal  indicators  in  order  to  decide  which 
process  should  be  done  first  because  the  two  sets  of  indi- 
cators are  unrelated.     — 

Before  a  design  is  performed,  growth  indicators  are 
determined  for  all  possible  growth  sites  and  the  largest 
indicator  is  selected.   If  no  transistor  grows,  removal  indi- 
cators are  determined  for  all  removable  transistors  and  the 
largest  indicator  is  selected.   In  most  cases,  only  one 
transistor  is  grown  or  removed  at  a  time  because  any  action 
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that  is  indicated  elsewhere  in  the  network  may  be  obviated 
by  the  action  on  the  transistor  with  the  largest  indicator. 
An  exception  to  this  procedure  is  when  a  differential  pair 
of  transistors  must  be  grown  together  to  retain  dc  balance. 

Our  first  approach  to  transistor  growth-removal'*"  con- 
sisted of  repeated  growth-removal  where  each  such  action  was 
followed  by  a  design  without  intermediate  supervision  by  the 
user.  Tlie  performance  function  remained  the  same  throughout 
the  process  and  the  weights  for  each  desired  iresponse  had  to 
be  chosen  very  carefully  in  order  to  avoid  overemphasis  of  a 
small  portion  of  the  desired  responses.  Many  runs  had  to  be 
made  with  different  weights  chosen  by  reducing  the  weight  on 
portions  of  the  desired  responses  that  were  receiving  the 
most  emphasis.  The  growth  sites  also  had  to  be  selected  for 
the  entire  run  and  their  best  locations  had  to  be  determined 
over  a  series  of  runs. 

It  is  difficult  to  specify  appropriate  sets  of  weights 
and  growth  sites  for  an  entire  run  so  a  new  approach  is  used 
that  allov/s  the  user  to  specify  new  sets  of  v;eights  and 
growth  sites  after  each  intermediate  design.   This  approach 
simulates  the  use  of  a  visual  display  and  an  interactive 
interface  v/hich  could  be  added  to  the  design  algorithm.   The 
flexibility  of  this  approach  allows  the  user  to  emphasize  a 
portion  of  the  desired  responses  for  one  or  more  growth- 
removal-design  steps  and  then  emphasize  another  portion  of 
the  desired  responses  at  a  later  step.   Thus  the  user's 
experience  may  reduce  the  effort  expended  by  the  computer. 
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The  examples  that  follow  were  run  with  intermediate  user 
interaction  on  the  computer  program  FROLIC. 

A  wideband  voltage  amplifier  is  the  basic  design  objec- 
tive of  the  first  example.   The  initial  configuration  shown 
in  Fig.  4-7 (a)  is  chosen  and  specifications  are  placed  on 
the  input  admittance,  output  inpedance  and  voltage  gain.   It 
is  convenient  to  express  these  specifications  in  terms  of  the 
hybrid  g-parameters  so  that  the  port  variables  of  the  amplifier 
are  given  by 

1^(0))  =  g^^(oo)l/j(w) 
and 

The  letter  A  in  Fig.  4-7 (a)  indicates  a  possible  transistor 
growth  site.   Gradients  for  voltage  gain,  A  ,  and  current  gain, 
A  ,  are  calculated  at  this  location.   If  a  gradient  component 
is  negative,  growth  is  indicated.   In  the  later  stages  of  the 
example,  several  growth  sites  will  be  available.   So  as  to 
prevent  unlimited  growth  or  removal,  only  one  transistor  is 
considered  at  a  time.   The  desired  responses  are  shown  in 
Fig.  4-7(b)-(d).   Performance  function  (3.6)  was  used  with 
p=2.   The  weights  were  specified  constant  over  frequency  at 
1  for  input  admittance  and  output  impedance  and  1000  for 
voltage  gain.   This  put  a  strong  emphasis  on  the  voltage  gain 
which  characterizes  a  voltage  amplifier.   Growth  indicators 
were  determined  and  a  common-base  stage  was  inserted.   After 
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Fig.  4-7 

(a)  Initial  configuration  for  voltage  amplifier  example 
(b)-{a)  Desired  (dashed)  and  actual  (solid)  responses 
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a  design,  the  resulting  network  and  responses  are  shov;n  in 
Fig.  4-8.   The  v/eights  were  then  adjusted  to  emphasize  the 
input  admittance  by  setting  the  weights  for  g    to  1000  and 
those  for  g    and  g    to  1.   Two  growth  sites  were  added  and 
growth  indicators  were  calculated.   A  common-collector  stage 
was  inserted  in  front  of  the  common-base  stage  and  a  design 
was  performed.   The  resxilting  network  and  responses  are  shown 
in  Fig.  4-9.   To  emphasize  the  output  impedance,  the  weights 
were  set  to  1000  for  g    and  1  for  g    and  ggi-   Two  new 
growth  sites  were  added  and  growth  indicators  were  calculated. 
A  common-collector  stage  was  grown  and  a  design  yielded  the 
network  and  responses  shown  in  Fig.  4-10.   The  weights  were 
adjusted  to  again  emphasize  the  voltage  gain  and  growth  indi- 
catoirs  were  calculated.   A  common-emitter  stage  was  grown  and 
a  design  yielded  the  network  and  responses  shown  in  Fig.  4-11. 
In  order  to  obtain  a  final  network  that  met  all  specifications 
more  equally,  all  weights  were  set  to  1  and  a  design  was  per- 
formed.  The  final  network  and  its  responses  are  shown  in 
Fig.  4-12.   The  voltage  gain  and  output  impedance  specifications 
are  met  satisfactorily  and  the  input  admittance  is  lower  than 
desired  which  is  quite  acceptable.   The  total  design  required 
35  seconds  of  IBM  360/65  computer  time. 

The  same  example  was  run  with  an  alternate  order  of  tran- 
sistor growths.   After  the  common-base  stage  v;as  grown,  empha- 
sis was  placed  on  output  impedance  instead  of  input  admittance 
as  before.   A  common-collector  stage  grew  behind  the  common- 
emitter  stage  and  a  design  yielded  the  network  and  responses 
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(a)  Conmon-base  stage  grew  to  satisfy  gain  requiremeni 
(b)"(d)  Desired  (dashed)  and  actual  "(solid)  responses 
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shown  in  Fig.  4-13.   More  growth  sites  were  added  and  input 
admittance  was  then  emphasized.   A  common-collector  was  grown 
and  a  design  yielded  the  network  and  responses  shown  in 
Fig.  4-14.   As  with  the  previous  example,  voltage  gain  v/as 
again  emphasized  in  a  run  and  then  equal  weights  v/ere  applied. 
The  final  network  and  responses  were  nearly  identical  to  those 
shown  in  Fig.  4-12  for  the  specifications  emphasized  in  the 
alternate  order. 

The  dual  of  the  voltage  amplifier,  a  current  gain  ampli- 
fier, is  now  obtained  using  the  same  interactive  procedure. 
Specifications  are  placed  on  input  impedance,  current  gain, 
and  output  admittance.   It  is  convenient  to  express  these 
specifications  in  terms  of  the  hybrid  h-parameters  so  that 
the  port  variables  of  the  amplifier  are  given  by 

l/j  (w)  =  hjj(co)I^(a)) 

and 

The  desired  responses,  h   (oo)  ,  h^.Cw),  and  h^^{ui)  ,    along  with 
the  actual  responses  and  initial  network  configuration  are 
shown  in  Fig.  4-15.   The  weights  of  input  impedance  and  out- 
put admittance  were  set  to  1  and  the  weights  of  current  gain 
were  set  to  1000  initially.   A  common-collector  stage  grew 
and  a  design  yielded  the  network  and  responses  in  Fig.  4-16. 
Input  impedance  was  emphasized  to  yield  the  network  and 
responses  shown  in  Fig.  4-17.   Emphasis  on  output  admittance 
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(a)  Common-collector  stage  grew  to  satisfy  gain  requirement 
(b)-(d)  Desired  (dashed)  and  actual  (solid)  responses 
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then  yielded  the  network  and  responses  shown  in  Fig.  4-18. 
The  current  gain  specification  was  then  emphasized  again  to 
yield  the  network  and  responses  in  Fig.  4-19.   To  obtain  a 
network  that  met  all  specifications  equally,  weights  of  one 
were  used  for  all  desired  responses  to  yield  the  network  and 
responses  of  Fig.  4-20.   The  complete  design  required  19 
seconds  of  IBM  360/65  computer  time. 

The  order  in  which  the  input  and  output  specifications 
were  emphasized  v/as  reversed  for  the  design  of  the  current 
amplifier.   The  network  and  responses  shown  in  Fig.  4-21  were 
obtained  when  output  admittance  was  emphasized  first.   After 
which  input  impedance  was  emphasized  to  obtain  the  network  and 
responses  shown  in  Fig.  4-22,   Further  emphasis  of  current 
gain  and  a  final  run  with  balanced  weights  yielded  the  same 
results  that  are  shown  in  Fig.  4-20. 

Since  transistor  removal  was  not  exercised  in  the  pre- 
ceding examples,  networks  were  constructed  that  required 
removal  to  demonstrate  the  removal  indicators.   A  fifth  tran- 
sistor was  inserted  into  a  version  of  the  four-transistor 
voltage  amplifier  and  a  design  was  performed  to  yield  the 
network  and  responses  shown  in  Fig.  4-23.   At  this  point,  all 
transistors  were  designated  removable  and  removal  indicators 
were  calculated.   Transistor  four  was  removed  and  replaced  by 
a  short  circuit  between  the  base  node  and  the  collector  node. 
A  design  was  performed  that  yielded  the  network  and  responses 
shown  in  P^ig.  4-24.   The  desired  responses  are  satisfactorily 
met  as  would  be  expected  since  the  transistor  that  was  removed 
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Fig.  4-19 

(a)  Conunon-emitter  stage  grew  to  satisfy  gain  requirement 
(b)-(d)  Desired  (dashed)  and  actual  (solid)  responses 
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Fig.  4-20 

(a)  Final  configuration  for  current  amplifier  example 
(b)-(d)  Desired  (dashed)  and  actual  (solid)  responses 
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(a)  Corniaon-base  stage  grew  to  satisfy 

output  admittance  requirement 
(b)-{d)  Desired  (dashed)  and  actual  (solid)  responses 
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Fig.  4-22 

(a)  Common-base  stage  grown  to  satisfy  input 

inpedance  requirement 
(b)-(d)  Desired  (dashed)  and  actual  (solid)  responses 
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Fig.  4-23 

(a)  Initial  configuration  of  voltage  amplifier 

with  extra  transistor 
(b)-(d)  Desired  (dashed)  and  actual  (solid)  responses 
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Fig.  4-24 

(a)  Final  configuration  of  voltage  amplifier 

after  transistor  had  been  removed" 
(b)-(d)  Desired  (dashed)  and  actual  (solid)  responses 
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was  the  one  that  had  been  inserted  and  it  was  replaced  prop- 
erly with  a  short-open  combination.   The  total  design  required 
6  seconds  of  IBM  360/65  computer  time. 

A  transistor  was  also  inserted  into  a  version  of  the  four- 
transistor  current  amplifier  previously  considered.   The  net- 
work, its  responses,  and  removal  indicators  for  each  transistor 
are  shown  in  Fig.  4-25.   Removal  of  transistor  two  is  suggested 
and  performed,  after  which  a  design  is  performed  to  yield  the 
netv/ork  and  responses  in  Fig.  4-26.   The  desired  responses  are 
again  m.et  satisfactorily.   The  total  design  required  4  seconds 
of  IBM  360/65  computer  time. 

The  growth  indicators  have  been  shown  to  accurately  pre- 
dict where  to  insert  transistors  and  in  v;hat  configuration  to 
put  them.   Although  the  principle  upon  v;hich  the  indicators 
are  based,  the  ability  of  a  transistor  to  supply  gain,  is 
simple,  it  is  used  here  to  solve  the  complex  problem  of  adding 
a  transistor  to  a  complicated  network. 

The  removal  of  a  transistor  from  a  network  is  not  as 
simple.   A  method  has  been  proposed  and  examples  have  been 
described  which  support  its  ability.   But  a  problem  does  exist 
with  the  method;  it  strongly  prefers  to  remove  common-emitter 
stages.   This  preference  is  not  understood  but  hopefully  an 
improved  method  will  be  available  in  the  future. 

For  a  complete  design  we  still  need  a  suitable  dc  bias 
netu'ork  to  supply  the  necessary  dc  collector  currents.   This 
design  typically  requires  the  ability  to  grow  additional 
transistors  to  meet  dc  specifications.   The  present  algorithm 


238 


B-lOO 

f  j=500   MHz       I  Ic'^-^  '^^ 
|l^I=.3   raA 

'^2 


3.3    Kfi 


fT=500    MHZ     l^cl  =  -3   ^^^ 


B=52  B=86 

g^3g  f j=492   MHz      f^.=500    MHz 

f.=500   MHz       |lcl  =  -9jJ^|lcl  =  -2    "^ 


T 


1   Mfi 


C 

T4 


^1 


i 


MQ. 


1   M.Q 


123    Kfi 


Tl 

T2 

T3 

T/| 

T5 

EB 

-5E9 

-4E10 

-lEll 

-5E10 

-7E9 

BC 

-4E10 

-7E13 

-2E12 

-3E11 

-6E6 

CE 

-4E8 

-2E12 

-5E10 

-7E9 

9E7 

(a) 


'11 

10  ^-J- 


10  * 


+ 


'21 
10^-'- 


10 


10 


-h 


10 


4-. 


3.. 


2-. 


+■ 


4- 


.005     .5      50 
Frequency  in  MHz 
(b) 


h 
10 


22 


10 


-5 


4- 


10 


-6 


-t- 


-t- 


.005     .5      50 
Frequency  in  MHz 
(c) 


4- 


,005    .5      50 
Frequency  in  MHz 
(d) 


Fig.  4-25 

(a)  Initial  configuration  of  current  amplifier 

with  extra  transistor 
(b)-(d)  Desired  (dashed)  and  actual  (solid)  responses 
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(a)  Final  configuration  of  current  araplifier 

after  transistor  had  been  removed 
(b)-(d)  Desired  (dashed)  and  actual  (solid)  responses 
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is  capable  only  of  adjusting  parameter  values  to  realize  the 
do  specifications.   Hence  the  circuit  designer  must  decide 
on  the  location  of  any  additional  transistors.   Further  in- 
vestigations will  hopefully  eliminate  the  necessity  of  this 
step.   These  amplifiers  can  be  realized  in  differential  con- 
figurations in  v/hich  the  dc  collector  currents  can  be  supplied 
by  simple  current  sources. 


CHAPTER  5 
SIMULTANEOUS  AC  AND  DC  DESIGN 

Although  separate  ac  and  dc  analyses  of  linear  networks 
are  efficient  and  informative,  separate  ac  and  dc  designs  are 
misleading  and  have  uncertain  value.   Some  network  elements 
influence  both  the  ac  and  dc  performance  of  a  network  so  that 
adjustment  of  such  elements  to  meet  ac  specifications  alone 
could  disturb  dc  perfor-mance .   The  ac  models  for  transistors 
and  diodes  are  dependent  upon  their  dc  operating  conditions. 
It  is  therefore  readily  seen  that  parameter  adjustment  to 
meet  ac  or  dc  specifications  alone  could  cause  deviation  from 
the  remaining  specifications. 

An  advantage  of  simultaneous  ac  and  dc  design  is  the 
ability  of  the  algorithm  to  reach  compromises  between  ac  and 
dc  specifications  if  necessary.   Some  specifications  require 
parameter  adjustments  that  would  keep  other  specifications 
from  being  met.   Quite  often  a  compromise  between  the  com- 
peting specifications  is  acceptable  or  even  desirable. 

Methods  have  been  discussed  in  the  previous  chapters  that 
were  selected  because  they  perform  ac  and  dc  design  efficiently 
and  accurately.   We  now  describe  how  these  methods  may  be 
combined  and  show  the  results  of  some  typical  designs. 
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Methods  Required  for  Simultaneous  AC  and  DC  Design 

Many  of  the  methods  used  to  perform  simultaneous  ac  and 
dc  design  are  identical  to  those  used  for  separate  ac  and  dc 
designs.   There  are  two  areas  that  require  revision;  these 
revisions  accomodate  the  design  algorithm  that  will  be  de- 
scribed later  in  this  chapter.   Specifications  are  included 
in  the  dc  design  that  temperature  desensitize  the  network 
lander  consideration.   In  order  to  implement  these  desensitizing 
specifications,  additions  are  made  to  the  performance  function. 
The  performance  function  and  the  parameter  gradients  become 
more  general  and  more  useful  due  to  the  modifications  that  are 
explained  in  this  section. 

Each  transistor's  operation  is  dependent  on  two  dc  varia- 
bles.  For  the  Ebers-Moll  model,  the  voltages  across  the  two 
diodes  are  considered,  and  for  the  hybrid-pi  model,  the 
collector  current  and  the  collector-to-emitter  voltage  are 
considered.   Each  diode's  operation  is  also  dependent  on  a  dc 
variable,  the  junction  voltage.   In  general,  these  dc  variables 
are  the  most  important  dc  variables  of  a  linear  amplifier.   If 
other  dc  variables  are  selected  as  representative  variables 
to  be  desensitized,  procedures  similar  to  those  that  follow 
can  be  developed.   In  the  discussion  that  follows,  a  general 
dc  variable,  Xj(r(T),T),  will  be  considered  as  it  is  repre- 
sentative of  all  possible  dc  variables. 

There  are  two  basic  approaches  to  desensitization  of  a 
variable  with  respect  to  temperature.   The  first  method  is  to 
desire  that  the  variable  assume  a  specific  value  at  all  tern- 
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peratures.   To  achieve  this  specification,  a  term  is  added 

to  the  pGrformance  function  for  each  dc  variable.   The  term  is 

n    T 
T    F 

E   =   ^   /    l/p{W,{T)  [X  ,(r{T)  ,T)/X,(T)-l]}PdT     (5.1) 
J=l   Tj        ^       ^  ~  ^ 

where  E^  is  the  contribution  to  the  performance  function  from 
temperature  desensitization, 
n   is  the  number  of  dc  variables  specified, 
Tj  is  the  initial  temperature  of  interest, 
T   is  the  final  temperature  of  interest, 
p  is  an  even  positive  integer, 
Wj(T)  is  a  nonnegative  weight  function, 
X  (r(T),T)  is  the  actual  value  of  the  dc  variable, 
r(T)  is  the  vector  of  network  parameters,  and 
X  (T)  is  the  desired  value  of  the  dc  variable. 

The  total  performance  function  can  be  written 

E  =  E.^+E^^+E^. 
AC   DC   T 

Although  this  approach  is  an  acceptable  method  of  de- 
sensitizing variables,  it  requires  the  user  to  select  a  value 
for  each  dc  variable  before  the  final  network  is  known.   Thus 
the  course  of  the  design  is  altered  to  achieve  the  specified 
values  of  the  dc  variables.   Even  if  this  method  is  useful 
for  some  design  situations,  an  alternate  method  is  required 
that  does  not  detract  from  the  other  objectives  of  the  network 
design. 

The  second  approach  to  temperature  desensitization  is  to 
select  a  standard  temperature  at  which  a  value  of  the  variable 
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will  be  determined  and  used  as  the  desired  value  for  the 
variable  at  other  temperatures.   Thus  a  specific  value  of  the 
do  variable  is  not  desired  and  the  desensitization  process 
can  be  accomplished  with  minor  influence  on  the  other  objec- 
tives of  the  design.   This  method  is  implemented  through  the 

performance  function  by  adding  a  term  of  the  form 
n   T 

^T  ==   ^   ^,   Vp{W^(T)[X^(r(T),T)/X^(r(T^),T^)-lJ}PdT      (5.2) 

where  T^  is  the  standard  temperature,  29  8°  Kelvin.   This 
approach  has  proven  to  be  very  satisfactory. 

In  Chapter  3  we  allowed  the  dc  bias  parameters  of  diodes 
and  transistors  as  well  as  the  nonbias  parameters  of  the  net- 
work to  vary  to  meet  the  ac  specifications.   For  simultaneous 
consideration  of  both  ac  and  dc  specifications,  the  dc  bias 
parameters  cannot  be  adjusted  arbitrarily  but  rather  assume 
values  which  are  dependent  on  conditions  existing  in  the 
network.   This  apparent  stumbling  block  can  be  removed  by 
making  a  simple  modification  to  the  ac  gradients  of  the  non- 
bias  parameters  which  incorporates  the  gradient  information 
for  the  bias  parameters. 

We  must  realize  that  many  of  the  nonbias  parameters  of 
the  network  directly  influence  the  values  of  the  bias  param- 
eters.  For  instance,  if  the  value  of  a  resistor  which  is 
connected  to  the  emitter  of  a  transistor  is  altered,  the 
value  of  the  emitter-base  junction  voltage  is  altered.   The 
interdependence  between  the  emitter  resistor  and  the  emitter- 
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base  junction  voltage  is  described  by  the  partial  derivative 

3V   /8R  . 
EB    E 

In  general,  a  bias  parameter  is  related  to  a  nonbias  param- 
eter by  the  partial-  derivative 

where  r   is  a  bias  parameter  and  r^,  is  a  nonbias  parameter. 

B  N 

Using  this  partial  derivative,  a  general  ac  gradient 
expression  can  be  determined  for  each  nonbias  parameter  that 
includes  gradient  information  for  itself  and  a  portion  of  the 
gradient  information  for  each  bias  parameter.   Thus,  the  ac 
gradient  component  for  each  nonbias  parameter  is  of  the  form 

3E   /3r   =  [dE      /dr    ]  +   T   [3E   /3r   ] [8r   /3r  ]       (5.3) 
AC    N        AC    N      j^j      AC    BI      BI     N 

I 

where  [3E   /3r  ]   is  the-  gradient  expression  heretofore  used 

to  directly  adjust  the  nonbias  parameter,  n   is  the  number  of 

bias  parameters,  and  8E   /r    is  the  gradient  expression 

AC   BI 

heretofore  used  to  directly  adjust  each  bias  parameter.   This 
expression  can  also  be  obtained  by  application  of  the  chain 
rule  of  differentiation  to  the  ac  performance  function, 

E   (r  ,r  (r  ))  . 
AC   N   B   N 

The  only  components  of  (5.3)  which  are  not  already  known 
are  the  partial  derivatives  of  the  bias  parameters  with  respect 
to  the  nonbias  parameters.   We  can  again  use  Tellegen's  theorem 
to  obtain  the  necessary  partial  derivatives  if  we  require  that 


146 


bias  currents  be  measured  through  voltage  sources  and  bias 
voltages  be  measured  across  current  sources.   If  we  refer 
back  to  (3.32),  rewritten  here  for  convenience 

n 

-Yai  \b   +yAV  (j)   =  -  y   G  Ar  , 
J   V^V  ^        II      j^j   J   J 

we  notice  that  the  partial  derivatives  are  easily  obtained 
after  an  analysis  of  the  original  network  and  n   analyses  of 
the  adjoint  network.   For  the  partial  of  a  voltage  bias  param- 
eter, v/e  set  all  voltage  and  current  sources  of  the  adjoint 
network  to  zero,  except  the  current  source  corresponding  to 
the  parameter.   Setting  the  remaining  current  source  to  unity 
yields   -  ' 

n 
AV^  =  Arg  =  -  I   GjArj 
J=l 

where  V^  is  the  bias  voltage  being  considered.   In  the  limit 
as  each  Ar   approaches  zero,  the  desired  partials  are 

3r„/8r,  =  -G,,    j  =  1,2,3, ...n.  (5.4) 

A  similar  result  can  be  obtained  for  a  current  bias  parameter 
if  all  adjoint  sources  are  set  to  zero  except  the  voltage 
source  corresponding  to  the  parameter  which  is  set  to  m.inus 
one. 

It  should  be  clear  that  the  dc  variables  used  to  temper- 
atiire  desensitize  a  network  are  identical  to  the  dc  bias 
parameters.   Although  this  correspondence  is  not  necessary, 
it  is  convenient  for  what  follows. 
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Since  we  have  added  a  temperature  desensitization  term 

to  the  performance  function,  gradient  components  for  that 

term  must  be  derived.   When  specific  values  are  desired  for 

each  dc  variable  as  shown  in  (5.1) ,  a  procedure  identical  to 

that  described  for  E    is  used.   When  the  value  of  the  var- 

DC 

iable  at  T   is  desired,  an  alternate  procedure  must  be  devel- 
0  ^ 

oped.   This  new  procedure  is  required  because  the  desired 
response  is  now  variable  also.   Previously,  all  desired  re- 
sponses were  constant  and  differentiation  of  the  performance 
function  yielded  an  expression  in  terras  of  partials  of  port 

voltages  and  currents.   Differentiation  of  (-5.2)  yields 

n^   Tp       • 

9E^/8r^  =  I      j         {Wj(T)  [Xj(r(T)  ,T)/X(r(TQ)  ,Tq)-1]}P~-^ 
J=l   Tj 

•W  (T){[9X  (r(T),T)/9r  ]/X(r(T  ) ,T  ) 

—        -   J  J  K  0     0 


-X(r(T)  ,T)/X(r  (T  )  ,T  )^8X(r(T  )  ,T  )/3r  }dT     (5.5) 


where 


8X(r  (T  )  ,T  )/9r   =  9r  /9r  . 
-  0  '  0     K     B'^   K 

The  only  quantity  of  (5.5)  that  needs  to  be  calculated  is  the 
partial  of  each  fixed  desired  response  with  respect  to  the 
network  parameter  of  interest.   These  partials  can  be  deter- 
mined in  a  procedure  identical  to  that  described  for  obtaining 
(5.4). 
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General  Algorithm  for  Simultaneous  AC  and  DC  Design 

The  basis  of  simultaneous  ac  and  do  design  is  to  deter- 
mine the  value  of  a  performance  function  that  reflects  the 
ac  and  dc  specifications  of  a  linear  amplifier  and  minimize 
it  by  means  of  an  optimization  algorithm.   As  indicated 
earlier,  the  optimization  algorithms  employed  require  deter- 
mination of  partial  derivatives  of  the  performance  function 
with  respect  to  the  designable  network  parameters.   The 
performance  function  and  parameter  gradients  are  determined 
by  the  appropriate  use  of  responses  of  the  original  network 
to  specified  excitations  and  responses  of  an  adjoint  network 
to  excitations  specified  above.   A  systematic  order  in  which 
to  perform  the  necessary  operations  will  nov;  be  discussed. 

Assume  a  linear  amplifier  is  to  be  designed  and  there 
are  both  ac  and  dc  specifications.   Transistors  are  used  to 
supply  the  amplifier's  gain.   Since  the  partial  derivatives 
of  the  bias  parameters  with  respect  to  the  nonbias  parameters 
are  required  for  both  the  evaluation  of  the  ac  gradients  and 
the  evaluation  of  the  dc  temperature  desensitization  portion 
of  the  performance  function,  these  partials  are  determined 
first.   A  single  dc  analysis  of  the  original  network  at  T 

0 

is  performed  and  the  values  of  all  of  the  dc  bias  parameters 
are  defined.   The  dc  adjoint  network  is  analyzed  n   times  and 
the  appropriate  responses  are  used  to  calculate  the  partial 
derivatives. 

We  choose  to  evaluate  the  ac  contributions  to  the  per- 
formance function  and  gradients  next,  even  though  the  dc 
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contributions  could  be  calculated  first.   For  the  ac  con- 
tributions, one  analysis  of  the  original  netv;ork  and  one 
analysis  of  the  adjoint  network  must  be  performed  at  each 
frequency.   The  contribution  to  the  performance  function  and 
excitations  for  the  adjoint  network  are  calculated  immedi- 
ately after  the  analysis  of  the  original  network.   The  adjoint 
analysis  is  then  performed  and  is  followed  by  the  gradient 
calculations.   Before  going  to  the  next  frequency,  the  partial 
integrals  are  augmented  with  the  appropriately  weighted 
contributions . 

Upon  completing  the  ac  calculations,  the  do  calculations 
are  performed.   For  each  temperature,  one  analysis  on  the 
original  network  and  n  +1  analyses  on  the  adjoint  network  are 
performed.   The  contribution  to  the  performance  function  from 
the  regular  specifications  and  the  corresponding  excitations 
to  adjoint  network  are  calculated  after  the  analysis  of  the 
original  network.   The  adjoint  network  is  then  analyzed  and 
the  gradient  contributions  related  to  the  regular  specifica- 
tions are  calculated.   For  each  do  bias  parameter,  an  addi- 
tional contribution  to  the  performance  function  is  determined 
and  a  corresponding  adjoint  analysis  is  performed.   Gradient 
contributions  related  to  each  bias  parameter  are  evaluated 
after  its  associated  adjoint  analysis.   Before  going  to  the 
next  temperature,  the  partial  integrals  are  augmented  with 
the  appropriately  weighted  contributions. 

The  performance  function  and  its  parameter  gradients 
have  now  been  evaluated  completely.   Control  is  passed  to  the 
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minimization  siibroutine  which  uses  the  gradient  information 
to  adjust  the  network  parameters.   Both  minimization  algo- 
rithms, conjugate  gradients  and  steepest  descent,  employ  a 
linear  search  follov/ed  by  cubic  interpolation.   This  approach 
requires  multiple  evaluations  of  the  performance  function  and 
its  parameter  gradients  for  each  linear  search.   After  deter- 
mining new  parameter  values,  the  minimization  subroutine 
returns  control  to  the  subroutines  that  evaluate  the  perfor- 
mance function  and  gradients .   The  procedure  again  starts  with 
the  evaluation  of  the  partials  of  the  bias  parameters  with 
respect  to  all  network  parameters.   The  order  of  operations 
is  identical  to  that  of  the  initial  calculations. 

The  examples  that  will  now  be  described  have  been  run  on 
FROLIC  which  is  depicted  through  flow  charts  in  the  appendix. 
A  simple  example  that  demonstrates  the  features  of  this  design 
algorithm  has  been  run.   The  desired  network  was  an  amplifier 
stage  biased  by  a  resistive  network.   The  initial  network  was 
an  RC  network,  shown  in  Fig.  5-1  (a),  with  a  centrally  located 
growth  site.   The  four  standard  bias  resistors  were  initially 
one  kilohm  and  the  remaining  resistors  were  ten  kilohms  each. 
Performance  functions  (3.3),  (3.6),  and  (5.2)  were  used  with 
p=2.   The  conjugate  gradients  minimization  algorithm  was  used. 
The  desired  responses,  shown  in  Fig.  5-l(b)-(c),  represent  a 
voltage  gain  of  100  and  minimum  dc  power  dissipation.   Growth 
indicators  were  determined  and  a  common-emitter  stage  was 
grown.   Upon  the  insertion  of  the  transistor  into  the  network, 
two  bias  parameter  temperature  desensitizing  specifications 


151 


12    V 


(a) 


'     0' 
10^  v-^ 


10  V 


1   V 


10' 


10' 


10< 


Frequency    in  Hz 
(b) 


10  mA" 


0  niA" 


293 


298 


303 


Temperature  in  °K 
(c) 


Fig.  5-1 


(a)  Initial  configuration  for  possible  amplifier  stage 

(b)  Output  ac  voltage  desired  (dashed) 
and  actual  (solid)  responses 

(c)  Do  supply  current  desired  (dashed) 
and  actual  (solid)  responses 
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were  invoked.   After  41  iterations  of  the  optimization  algo- 
ritl-na,  158  function  evaluations  had  been  performed  and  35 
seconds  of  IBM  360/65  computer  time  consumed.   The  resulting 
network  and  its  responses  are  shown  in  Fig.  5-2.   The  param- 
eters and  constants  of  the  final  network  are  shown  in  Table 
5-1.   The  desired  specifications  had  been  met  satisfactorily. 

A  more  complicated  example  involving  a  differential  pair 
of  transistors  and  a  transistor  current  source  was  also  designed. 
The  initial  netv7ork  and  its  desired  responses  are  shown  in  Fig. 
5-3.   Till  resistors  were  initially  one  kilohm.   The  desired 
responses,  performance  function,  and  minimization  algorithm 
were  basically  identical  to  those  used  in  the  previous  exairiple. 
The  differential  pair  v/as  grown  upon  calculation  and  investi- 
gation of  the  growth  indicators.   In  this  example,  tv.-o  tran- 
sistors were  allowed  to  grow  at  once  as  the  network  required 
a  balanced  pair.   Since  the  network  contained  four  transistors, 
eight  temperature  desensitizing  specifications  were  invoked. 
After  .15  iterations,  32  function  evaluations  had  been  performed 
and  17  seconds  of  IBM  360/65  computer  time  had  been  consumed. 
The  resulting  network  and  responses  are  shown  in  Figs.  5-4 
and  5-5.   The  parameters  and  constants  of  the  final  network 
are  shov/n  in  Table  5-2.   Again  the  desired  results  were  obtained, 

An  eleven-transistor  video  amplifier,  the  Fairchild 
ViA.75lC,  is  a  good  test  for  the  design  algorithm.   Since  this 
network  incorporates  feedback,  the  automatic  frequency  selec- 
tion algorithm  was  used  to  guarantee  that  all  Integrals  were 
evaluated  accurately.   Transistor  growth  was  not  exercised 
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(a)  Final  configuration  of  one-stage  transistor  amplifier 
(b)-(e)  Desired  (dashed)  and  actual  (solid)  responses 
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Table  5-1 

Parameters  and  Constants  of  Final  Configuration 
of  One-Stage  Amplifier 
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(a)  Initial  configuration  for  differential  amplifier  example 
(b)-(d)  Desired  (dashed)  and  actual  (solid)  responses 
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for  temperature  desensitization  specifications 
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Table  5-2 

Parameters  and  Constants  of  Final  Configuration 
of  Differential  Amplifier 
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in  this  example.   The  desi.red  responses  represent  an  ac  volt- 
age gain  of  800  and  a  minimum  do  power  dissipation.   All 
resistors  were  started  at  one  kilohm  and  transistor  betas 
and  f  s  were  started  at  50  and  500  MHz,  respectively.   Asso- 
ciated with  the  11  transistors  were  22  bias  parameters,  each 
of  v/hich  was  used  as  a  temperature  desensitizing  specification, 
Performance  functions  (3.3),  (3.6),  and  (5.2)  v/ith  p=2  and 
the  conjugate  gradients  minimization  algorithm  were  used. 
After  44  linear  search  iterations,  117  function  evaluations 
had  been  performed  and  1,018  seconds  of  IBM  360/65  computer 
time  had  been  consumed.   The  final  network  and  its  responses 
are  shown  in  Figs.  5-6  and  5-7.   The  parameters  and  constants 
of  the  final  network  are  shown  in  Table  5-3.   The  resultant 
network's  resistors  are  larger  than  those  of  the  ijA751C  be- 
cause of  the  minimum  dc  current  'specif ication.   If  a  parasitic 
capacitance  were  added  to  each  resistor  and  its  value  v;as 
proportional  to  the  resistance,  the  resistor  values  V70uld  have 
been  smaller  and  closer  to  those  of  the  viA75lC.- 

A  design  procedure  has  been  proposed  and  demonstrated 
that  simultaneously  achieves  ac  and  dc  specifications  on  a 
linear  amplifier.   The  procedure  has  been  implemented  in  a 
computer  program,  FROLIC,  which  incorporates  techniques  that 
are  efficient,  accurate,  and  versatile. 
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CHAPTER  6 
CONCLUDING  REMARKS  AND  SUGGESTIONS  FOR  FURTHER  RESEARCH 

An  automated  network  design  algorithm  has  been  described 
which  simultaneously  considers  ac  and  dc  specifications. 
Advanced  techniques  have  been  assembled  to  best  accomplish 
the  task  of  designing  the  large,  complex,  integrated  circuits 
in  use  today.   Examples  have  demonstrated  the  ability  of  the 
algorithm  to  efficiently  achieve  desired  specifications. 
Some  additional  work  which  v/ould  enable  the  algorithm  to  be 
more  flexible,  miore  accurate,  and  more  efficient  is  described 
below. 

Improved  models  for  the  physical  circuit  elements  which 
account  for  parasitics  are  necessary.   In  particular,  resistor 
models  need  to  include  capacitive  and  inductive  effects. 
Bipolar  transistor  models  require  a  collector  current  depend- 
ent  beta. 

Consideration  of  additional  types  of  circuit  elements 
would  broaden  the  scope  of  the  algorithm.   Elements  such  as 
field-effect  transistors  and  RC  structures  are  often  used  in 
integrated  circuits.   Attention  should  also  be  given  to  often- 
used  combinations  of  elements  such  as  current  sources  and 
differential  amplifiers.   As  new  elements  and  groups  of  ele- 
ments are  considered,  models  and  gradient  expressions  for 
these  components  must  be  determined. 
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Growth  and  removal  of  bipolar  transistors  to  meet  ac 
specifications  have  been  investigated.   Further  work  should 
be  done  in  this  area  with  consideration  given  -to  the  do 
characteristics  of  transistors.   Growth  and  removal  in- 
dicators should  be  determined  that  consider  the  total  per- 
formance of  transistors.   Possible  growth  and  removal  of 
diodes  and  other  elements  that  may  be  added  should  be 

explored. 

Tolerances  related  to  the  various  steps  of  the  reali- 
zation of  an  integrated  circuit  make  it  difficult  to  obtain 
the  exact  network  that  was  specified  by  the  design  algorithm. 
These  tolerances  should  be  accounted  for  during  the  design 
process  so  that  resulting  networks  are  capable  of  meeting  the 
desired  specifications  without  obtaining  the  exact  values 
designated  for  all  of  the  parameters. 

Temperature  desensitized  ac  gain,  a  desirable  quality  of 
an  amplifier,  should  be  incorporated  in  the  design  specifi- 
cations.  Complete  consideration  of  this  specification  would 
require  ac  analysis  of  the  amplifier  at  a  set  of  frequencies 
which  represents  the  frequency  range  for  each  temperature  of 
a  set  that  represents  the  temperature  range.   This  time- 
consuming  process  would  have  to  be  carried  out  for  each 
evaluation  of  the  ac  contributions  to  the  performance  function 
and  gradients.   A  condensed  method  which  may  closely  approxi- 
mate the  above  process  is  to  analyze  the  amplifier  at  a  set 
of  frequencies  for  a  standard  temperature  and  then  analyze  at 
a  set  of  temperatures  for  a  standard  frequency. 
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As  new  options  are  added  to  the  automated  network  design 
algorithm,  the  ability  to  produce  nev?  and  improved  circuits 
is  increased.   The  rate  at  v/hich  these  circuits  are  designed 
can  be  increased  significantly  if  results  of  previous  designs 
are  available  to  the  algorithm.   A  good  initial  network  is 
important  in  many  cases  and  a  helpful  hint  is  often  needed 
in  the  midst  of  the  design  process.   This  feature  can  be 
realized  in  two  ways:   devote  large  amounts  of  computer 
memory  to  storing  all  significant  circuit  configurations  or 
allov.T  interaction  by  a  human  circuit  designer.   Since  large 
amounts  of  time  would  be  required  to  maintain  and  search  a 
storage  area  full  of  circuit  configurations ,  the  more  prac- 
tical approach  is  to  develop  an  interactive  system  that  allows 
a  designer  to  monitor  the  progress  of  a  circuit  design. 
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